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ABSTRACT 


The present paper is a reflexion 
on the underlying criteria that 
preside to the selection of che- 
mical elements for specific 
functions in living organisms, 
and a brief review of the physi- 
cal-chemical principles on 
which their uptake from the 
environment, transport and 
trapping is based, As shown, 
not only thermodynamic but 
kinetic and feed-back control 
conditions must be fulfilled, so 
that life chemistry can not be 
regarded simply as a multiple 
set of reactions tending to a 
state of equilibrium with the 
environment (which would 
mean death), but as a dynamic 
process of continuous inter- 
actions of living species with 
the environment, whereby the 
environment is changed and 
the changes are reflected back 
in the chemistry of the living 
species themselves. 


SUMÁRIO 


O presente artigo constitui uma 
reflexão sobre os critérios sub- 
jacentes à selecção dos ele- 
mentos químicos para funções 
especificas nos organismos 
vivos, € uma revisão breve dos 
princípios fisico-químicos nos 
quais se baseia a sua captura 
do ambiente, transporte e fixa- 
ção interna. Tal como se des- 
creve, há que satisfazer condi- 
ções termodinâmicas, cinéticas 
e de controle retroactivo, e por- 
tanto a vida não pode ser 
olhada simplesmente como um 
conjunto de reacções tendendo 
para um estado estável de 
equilibrio com o ambiente (o 
que equivaleria à morte), mas 


THE SELECTION OF 
CHEMICAL ELEMENTS BY 
BIOLOGICAL SYSTEMS -— 

SOME THERMODYNAMIC AND 
KINETIC ASPECTS 


J.J.R. Fraústo da Silva 


Centro de Quimica Estrutural, Instituto Superior Tecnico 


1. Biological elements and their functions 


The majority of present day living organisms require rather strict environmental 
conditions: temperature from O to about 40ºC, pressure of the order of 1 atm, 
salinity up to about 4%, pH in the range 4 to 9, redox potential in the region of 
stability of water (-0.4V to +0.8V at pH 7, vs. the standard H*/H, electrode), 
water activity between 0.7 and 1, etc. A few exceptions, especially the so-cal- 
led extremophiles, confirm these general conditions — see Table |. 


Table | 
Minimum and maximum values of environmental factors of present day life 


Lower limit 
-18"C (fungi, bacteria) 


Factor 
Temperature 


Upper limit 


+110ºC (sulphate reducing bacteria 
under 1000 atm hydrostatic pressure) 


Redox potential -450 mV at pH 9.5 (sulphate +850 mV at pH 3 (iron bacteria) 
reducing bacteria) 
pH O (Acontium velatum, fungus D, 13 (Plectonema nostocorum) 


Thiobacilus thigoxidans) 


Hydrostatic pressure Essentially zero 1400 atm (deep-sea bacteria) 


Salinity Double distilled water Saturated brines (Dunalielia, 
(heterothropic bacteria) halophyllic bacteria, etc.) 
Water activity 0.65-0.70 (Aspergillus glaucus) Essentially 1 


Source: J.R. Valentyne in G. Mamikunian and M.H. Brigs (eds), “Current aspects of exobio- 
logy, Pergamon Press, London 1965, pp. 1-12. See also M.T. Madigan and B.L. Marrs: 
Scientific American, April 1997, pp. 66-71. 


The required environmental conditions have their counterpart in the neces- 
sity of an equally strict internal regulation of processes, and this means time 
and space organisation of the biological molecules and of their reactions. 
Such is the scope of the chemistry of life but, differently from the classical 
tradition, we regard it as the interplay of not only a set of small and of very 
large organic molecules, but of these and a group of about 25 chemical ele- 
ments, free as ions or combined, organised in the biological space and time. 
These elements, considered as “essential” for animal or vegetal life, are 
identified in Figure 1, the Periodic Table of the Elements, and a first conclu- 
sion is immediate: that the elements essential for life are the lighter ele- 
ments, metals and non-metals of atomic number less than 35. Amongst 


required elements only molybdenum, tungsten, iodine and tin have higher 
atomic numbers and none of these seems to be generally required by all 
living species. 
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Fig. 1 The essential elements in life. 


On the other hand, the biological elements seem to have been selected from 
practically all groups and sub-groups of the Periodic Table (the only exceptions 
are groups 3 and 4, besides that of the inert gases) and this means that practi- 
cally all kinds of element properties are associated with life processes, within 
the limits imposed by external constraints. The reason for the selection of 
these elements appears to have been determined by abundance in the earth's 
crust and availability (solubility) in water under primitive anaerobic conditions. 


Before showing how can these element properties be put to use and eventu- 
ally adapted and perfected by evolution, we must recall that oxygen and 
hydrogen are, to a large extent, combined with each other in the form of 
water, and that most life chemistry takes place in aqueous media. This fact 
implies that, in general, only processes compatible with the presence of this 
solvent are possible, restricting the utilisable range of redox potentials, pH 
and temperature. 


Bulk metals or alloys, for instance, cannot be used by biological systems 
since they are not stable in water, and processes depending on metallic 
electronic conduction which might be valuable, e.g. to send messages, must 
be carried out by other chemical agents or systems. 


The same can be said for temperature constraints, hence the high energy 
requiring operations of industry used to obtain products that biological 
systems also need, e.g. NH, from Ns, must be carried at normal physiologi- 
cal temperatures and this is only possible by using extremely sophysticated 
catalytical processes developped over thousands of millions of years. 


As to the functions that the various chemical elements can play in biology, 
depending on their properties eventually modulated and adjusted by specific 
combinations or situations, the following list of well known examples gives an 
idea of the extension of the involvement: 


a. Biological molecules are covalent compounds of the bulk elements (car- 
bon, hydrogen, oxygen and nitrogen), sometimes with sulphur or phos- 
phate groups and in a few cases with some minor elements, e.g. sele- 
nium (in glutathione peroxidase) or iodine (in thyroid hormones). 

b. The hard structures of biology — bone, teeth, shells, etc. — are largely 
mineral, containing silica or calcium salts. 


c. The release of oxygen from water, that maintains our present atmos- 
phere, is absolutely dependent on manganese (component of an O, — 
evolving enzyme). 


d. Oxidative metabolism is dependent on iron and copper catalysts. 


como um processo dinâmico € 
continuo de interacção dos 
organismos vivos com o ambi- 
ente, através do qual este sofre 
alterações e as alterações se 
repercutem na quimica interna 
dos próprios organismos vivos. 


e. The capture of nitrogen from nitrate or dinitrogen relies totally on molyb- 
denum containing catalyst centres. 


f. The synthesis of DNA and RNA and the activation of water in condensa- 
tion and hydrolysis reactions, e.g. in digestion, depend on the intervention 
of more than 200 zinc enzymes (DNA synthesis also requires either 
cobalt or iron in ribonucleotide reductase). 


g. Phosphate transfer, essential in bioenergetics, is almost invariably linked 
with magnesium activation. 

h. The capture of light in photosynthesis also depends on a magnesium pig- 
ment, chlorophyll. 

|. Electrolytic connection in cells and between cells, e.g. nerve messages, 
rests largely on the movements of sodium, potassium and calcium. 


|. Mechanical changes in cells are initiated by calcium fluxes and the same 
is true for triggering in digestion, blood-clotting, fertilization, muscle con- 
traction, cell-surface interactions, etc. 


|. Transport of oxygen requires iron or copper proteins, etc. etc.... 


These are just a few illustrative cases and many more could be presented, 
but they are enough to show that a wide range of functions is covered — 
structural and mechanical functions, electrochemical control, acid-base and 
redox catalysis, generation, transmission and storage of energy, transmis- 
sion of messages and triggering effects, reproduction, etc. 


Furthermore, one can identify some specialisation in these various functions: 
iron and copper are associated with oxidation and electron-transfer, zinc with 
the metabolism of water, molybdenum with the metabolism of nitrogen, cal- 
cium with triggering and other mechanical effects, sodium and potassium 
with electrolytic conduction, etc. 


Such specialised roles are related to the properties of the corresponding ele- 
ments, but there are also special reasons for some individual preferences. 
Copper, for instance, could, in principle, replace zinc in hydrolytic reactions 
since it is as good a Lewis-acid as zinc is, but there may be some risks in its 
utilisation due to its easy variation of oxidation state and to the eventual for- 
mation of free-radicals. 


Although the list of examples given stresses the role of the metals, several 
non-metals are also essential for life, besides those involved in the organic 
biological molecules, i.e. C, H, N, O, S and phosphate groups. That is the 
case of silicon, with structural functions; selenium, involved in some redox 
reactions; chlorine (as Cl), involved in electrolytic conduction; iodine, cova- 
lently bound in thyroid hormones; and boron, with a limited role in certain 
higher plants. Still, it may be considered that the non-metals are more in the 
realm of classical biochemistry, even if this distinction is rather artificial, and 
in the present paper we will concentrate mostly on metal ions, particularly on 
the processes used by biological systems to select some that are assential 
and reject others that are unnecessary or toxic. 


2. THE UPTAKE AND PARTITION OF THE BIOLOGICAL ELEMENTS 


We shall now examine how biological systems can take up and distribute the 
elements they need for their metabolic functions. This will require us to dis- 
cuss in somewhat greater detail the mechanisms used to achieve selectivity 
and the physico-chemical basis of the competition of the elements for biolo- 
gical molecular or solid structures. 


We will be concerned with the mechanisms operating at cell level and leave 
aside the details related to the environmental conditions. For all that matters, 


NH 


we will consider the elements in aqueous solution, in a state determined by 
the pH and the redox potential (high). There might be certain ligands L, pre- 
sent and the elements M considered may be party complexed as shown in 
Figure 2. 


OUTSIDE CELL “ INSIDE CELL 


| Cell emma 


ÉS 
L - q— 1 
do 
/ W ai neto Trap 
| ML, ML, E 
o M Qu 
kd M L, 
E: high E z EE RS o E* low 
pH variable pH=7.0 
CELL MEMBRANE | 


Fig. 2 Ligands L,,L,, and L, are synthesized by the cell. Ligand L, represents extracellular 
complexes, 


This Figure shows, in general terms, the parameters that must be taken into 
account in both uptake and control of intracellular concentrations of trace 
metais. 


The free external concentration of M is determined by abundance, availability 
and the environmental factors pH, E, and the stability constants for any ML, 
complexes. The internal concentration is determined by: 


a. The complexation by selective carrier ligands L.. 


b. The capture by ligands L, at the cell surface that can remove M from 
solution. 


c. Diffusion through selective channels due to concentration gradients esta- 
blished accross the external membrane of the cell by passive, assisted or 
active (energy-coupled) processes. 


d. Binding to internal polymers from which exchange of metal is very slow. 


Hence, uptake is controled by metabolic activity, that is responsible for the 
synthesis of L,, L, and L, and for energy-coupled processes requiring in- 
vestment of cellular-energy, which are all time-dependent phenomena. 
Clearly, one must consider not only equilibrium thermodynamic parameters 
related to metal ligand-binding, but rates of input from the environment, rates 
of transport accross cell membranes, rates for inclusion in cellular traps L, 
and even rates of eventual exit from the cell, The selection of elements (and 
their rejection) are therefore complex processes whose mechanisms involve 
not only equilibrium thermodynamics but also reaction kinetics, which must 
be considered together. 


With the exception of studies on the transport of iron into microbial cells and 
on the transport of alkali and alkaline earth-metals such as Na*, K* and Ca”, 
there in still very little basic understanding of the processes involved in the 
cellular uptake of essential metal ions, such as Cu**, Ni”, or Cr*. This area 
has been rather neglected in biochemistry and, consequently, we do not 
understand well several important problems of present concern, e.g. anoma- 
lous reasons for nutrition deficiencies or the regulation of concentrations of 
transition metals below toxic levels in polution situations. Our theories have 
been built on the assumption that the intracellular compartments are in equi- 
librium with the cell environment and until recently we haven't realized that 
this is a typical example of a non-equilibrium process. 


Before considering the components of this process it is useful to explore fur- 
ther the most common device that biological systems use to divide space 
and separate different phases: the biological membranes. 


Indeed, biological membranes are outstandingly important but not very well 
known systems: not only they divide space but they determine mobility in 
and accross themselves and can incorporate a number of energy using 
systems that can move elements from one compartment to another either as 
simple ions, small chelates, protein complexes or even crystallites. On the 
other hand, membranes are tri-dimensional assymetric bodies, including in 
this assymetry the distribution of their component systems. 


The energy using systems are proteins which act as pumps fixed in the 
membrane or whereby the membrane can invaginate into the cell to form a 
new vesicle or coalesce with an internal membrane to discharge the contents 
into a existing compartment. Hence they not only separate compartments but 
they can also regulate the trafic of biological substances — ions, ligands or 
complexes — between compartments. 


These compartments may have different conditions, namely pH and redox 
potential, hence the situation of those biological substances may be changed 
in the process, i.e. they will be energised relatively to a certain reference 
state. For instance, the state of a ligand carrying two thiol groups can be 
(-SH)., (-S), or -S-S- depending on the values of pH and redox potential in 
the compartment into which the ligand is pumped or in which is synthesised 
in variable amounts by cell metabolism, as an example, free thiol groups are 
not stable outside cells since the redox potential is high. Similarly, some 
metal complexes dissociate in compartments where pH is too low. 

These observations point to the logical deduction that the distribution and 
localisation of the various elements in biological systems is necessarily ass- 
ymetric, even at the cell level. We will come to this aspect later. 

Before proceeding it is useful to stress that the access of charged species to 
cells is limited by their hydrophyllic character since membranes are hydro- 
phobic. Neutral molecules, e.g. HO, CO,, No, NH,, O,, B(OH), or Si(OH), 
diffuse freely through membranes and the problem here is to prevent back 
diffusion. Obviously, these molecules or their elements must be trapped 
inside cells by various means, usually by redox-reactions followed by con- 
densation of the products with internal molecules, forming thermodynamically 
labile but kinetically stable bonds (covalent). 


lonic species, e.g. metal ions, must be dealt with by different process which 
must depend on the individual properties of such species. 

The first aspect to consider is the real binding of these species to various 
ligands, which is related to the binding constant K,, IN any one compart- 
ment. Since the process is open to competition and most biological ligands 
are involved in acid-base equilibrium, the constants to be considered are the 
so-called “effective” or “conditional” stability constants. Even so, these effec- 
tive constants are just a first approximation since the biological ligands of 
interest are in many cases rather complicated polymers, e.g. proteins, for 
which the necessary metal-complex stability constants are not known in real 
situations. The values used for K,q, are for single “model-compounds”, usu- 
ally small organic molecules containing an equivalent set of donor atoms and 
desirably, as much as possible, in a geometry equivalent to that of the site of 
the natural ligands. This site is usually distorted and the energetics of the 
real situation are different from those of the “model”, hence the K,q, Values 
are sometimes rather rough approximations. On the other hand, the proteins 
are large and flexible molecules that fold in different ways on complexing a 
metal ion or a complex, originating several new internal interactions among 
which those of the metal ion with the second and third coordination spheres. 
This effect may originate considerable differences in stability reaching seve- 
ral orders of magnitude, see later, With these restrictions in mind we can 
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now examine the various factors that affect the selective uptake of elements by 
the biological systems. 


3. UPTAKE OF METAL IONS AT EQUILIBRIUM 


The examination of the different properties on which selective uptake can be 
based may start with the inspection of the tendencies and trends exhibited by 
metal ions as regards the formation of complexes, influencing the values of 
Kay): Some general characteristics are shown in Tables Il to IV. 


Table Il 
Some general characteristics of metal ions present in biological systems 


Characteristics Na,'K* Mg?” ,Ca* Zn” Ni” Fe,Cu,Co,Mo,Mn 
Stability of the Very low Low or medium High High, but not so 
much in the cases 
complexes of Mn* and Fe” 
Preferred donor Oxygen Oxygen Sulphur (or nitrogen) Nitrogen of sulphur 
atoms (and oxygen in the 
case of Fe” and 
Mn*) 
Predominant +1 +2 +2 More than one state 
oxidation state 
Mobility in biological Very mobile Semi-mobile Static Static (semi-mobile 
media in the cases of Mn” 
and Fe”) 
Table Ill 
Ligands preferred by different metal ions in simple coordination compounds 
Metal ion Ligands 
Na K' Oxygenated ligands, neutral or of low charge (-1) 
Mg” Carboxylate, phosphate and polyphosphate 
Ca” Carboxylate, phosphate (less than Mg”), neutral O-donors 
Mn” Carboxylate, phosphate (nitrogen) 
Fe” Thiolate (e.g. cysteine), unsat. amines (particularly porphyrines) 
Cr”, Mn” Phenolate, (e.g. tyrosine), hydroxamate, hydroxide, 
Fe” Co” carboxylate, porphyrins, polyrroles 
Cu Ni” Thiolate (e.g. cysteine), unsat amines, polypyrroles 
Cu” Amines, imidazole, ionised peptide, thiolate 
Zn” Amines, imidazole, thiolate 
Table IV 
Preferred geometries in simple coordination compounds (*) 
Metal ion Preferred geometries 
Cu” Mn” Tetragonal > 5 — coordination > tetrahedral 
Ni” Octahedral > others 
Co* Octahedral > tetrahedral > others 
Znº Tetrahedral > octahedral > 5 - coordination 
Mn Octahedral > others 
Fe” Co*,Cr* Octahedral > others 


() N.B. all metal ions are in high-spin states and liganging atoms are smal O,N donors. 


On the other hand, the stability of the complexes formed by the different types 
of ions follows sequences which are, to a large extent, independent of the 
ligand. 


For the common biological monovalent ions the binding constants follow the 
order Cut > Na*, K* and for the divalent ions one has the so called Irving — 
Williams or “natural” stability order: (Mg”*), Ca” < Mn” < Fe” < Co” <« Ni? < 
Cu? > Zn?. Few stability constants are available for trivalent ions, but, in 
general, we expect and find the order Co(lll) (low-spin) > Mn(Ill) > Fe(lll) > 
Cr(ll) > AI(II). 

Naturally, trivalent ions usually form stronger complexes than divalent ions 
and these form stronger complexes than monovalent ions. Sodium and 
potassium, particularly, form relatively weak complexes, hence they 
exchange and diffuse rapidly as simple tons and can only be retained effecti- 
vely by physical barriers to ionic diffusion, i.e. by membranes between 
aqueous domains into which these ions must be “pumped” using metabolic 
energy. For magnesium and calcium, as well as for the transition metals and 
Groups 11 and 12 it becomes increasingly possible to retain them by binding 
in polymers (proteins, polysacharides), but their final site often is required to 
be in a state ffom which exchange is very slow, 1.e. for which the thermody- 
namic binding strength and/or the kinetic activation energy for escape is 
made very large, in order that they can be of use as fixed reaction centres. 
These ions represent a gradation from Group 1 to the non-metals or from 
electrostatic to covalent binding. 


In these general comments we have already stressed two selection factors — 
charge and ionic size — but others can be equally important, as it is obvious 
from Tables Il to IV: preference for certain donors, preference for a particular 
coordination symmetry, stabilisation by spin-pairing. Other possibilities 
include cooperative binding of clusters and changes on going from aqueous 
solution to the inside of proteins (transfer coefficients). 

Naturally, the amount of metal bound to a ligand is proportional to its con- 
centration in the medium ([ML]..K,y .[M]), hence relative concentrations of 
metals and ligands is also a factor that determines uptake. 


We will comment briefly on these vanous factors to show how they are put to use: 


a. Charge — Charged ligands of smaller size prefer the highly charged class 
“a” cations and this means that appropriate donors for tri-positive ions are 
small first row anions such as >N', RO' and F. Since F' cannot be combi- 
ned in an organic group and the concentration of >N' is small at pH 7, this 
leaves RO' (phenolates, hydroxamates, enolates) as the adequate donor 
for tri-positive metal ions, e.g. Fe”, Mnº* and probably Cr*. This is in fact 
true. 

b. Size — Since the order of stability for electrostatic binding in that of the 
ionic potencial, one expects the following orders: 


Li'>Na'>K'>Rb'>Cs* 
Be**>Mg''>Ca“>Sr'>Ba” 


Due to stereochemical constraints to the complexation of small ions by 
polydentate ligands, these orders are seldom found. The inversions can 
only be explained in terms of energy balances involving the hydration 
energies of the metals and the ligands and the thermodynamic effects in 
the complexation reactions, hence no typical pattern is obvious. 


Typical examples of selection by size are those of the bacterial antibiotics 
nonactin and valinomycin which are used to preferentially capture K*, and 
of chlorophyll which takes up Mg” selecting against Ca? (although other 
factors may be involved). We also find many calcium-binding proteins, 
e.g. calmodulin, parvalbumin, thermolysin, etc. etc, with specific sites for 
Ca“ (Mg2+ not bound) and some Mg*-activated enzymes, e.g. phos- 
phomutases, with coordinating sites for Mg? in which Ca? binds poorly 
(although in this case it is not only a matter of size). 


[oi 


im 


The best examples of selection of metal ions based on ionic size alone are 
probably those of the specific protein channels for Na*, K* and Ca**, but 
here no binding is necessary and is indeed undesirable. The mechanism is 
by diffusion or pumping and it is Kinetical rather than thermodynamical. 


Type of donor — As it is well known, A-group metal ions prefer O-donor 
ligands whereas B-group metal ions prefer to coordinate to sulphur or 
nitrogen. Furthermore, better donors, which form more covalent bonds, 
intensify the differences of stability between the complexes of the 
B-group of metal ions, as shown in Figure 3. This Figure, which illustra- 
tes the Irving-Williams' order of stability, shows that the divalent transition 
metals form complexes that may differ in their thermodynamic stability by 
several orders of magnitude if the ligand atoms are appropriately chosen. 
This is why metal ions such as Cu*, Ni”, Cd” are so dangerously toxic 
when in excess in the biological environment; they will be preferentially 
complexed rather than other essential metals, e.g. Fe” or Zn”. 


| 
TO + 
| NH>CHsCH;S) 
e 
| | 
Ea 
m 6| A 
o A 
4 
OM 
x 
il -NHsCHs5CO, 
O 
l Í ain. e JR RD a PR A J 


Ba A Ca Mg Mn Fe Co Ni Cu Zn 


12 Fig. 3 The variation of stability constants, K, for the complexes of M** ions of the Irving- 
-Wiliams series, ox, oxalate; en, ethylenediamine. 


a. 


Preference for coodination symmetry — Itis a well known fact of coordina- 
tion chemistry that metal ions have stereochemical coordination preferen- 
ces. For instance Cu” prefers tetragonal coordination arising from the 
Jahn-Teller effect; on the contrary, Znº* has a demand for tetrahedral 
geometries found in most biological zinc sites, e.g. in carbonic anhydrase 
and alcohol dehydrogenase. Another case of different stereochemical 
preferences is that of the pair Mnº*t and Feº*; the first ion prefers tetrago- 
nal sites and the second prefers tetrahedral sites. 


By combining stereochemical preferences and types of ligands we can 
build a Table of selectivity in which the favoured sites become obvious — 
Table V — and several examples of these are found in biological mole- 
cules such as the blue-green Cu**/Zn**, the red Mnº* and the yellow Fe* 
superoxide dismutases. 


Table V 
Table of selectivity for Cu”, Zn”, Mn and Fe (selection by adequate 
combinatiion of donor atom and geometry of site) 


Geometry of site Donor atom 


Phenolate RO lonised peptide>N' Thiolate RS Nitrogen N 


Tetrahedral Fe” ” Zn” Fe” Zn 
Tetragonal Mn* Cu E: Cuê 


e. Spin-pairing stabilisation — For transition metal ions there is a further pos- 
sibility to increase selectivity based on the free-energy change caused by 
spin-pairing in higher ligand-fields. The ions Coº* and Fe” (both dº type) 
are favoured and we expect and find that strong N/S donors bind 
Co**>Mnº%>Fe*>Cr** and Fe**>Co?*>Mn*. 


The selection of Fe”? preferably to Cu** in cytochromes and of Co** for 
vitamin B,, in the presence of the other transition metal ions may be, at 
least partly, based on this effect. 


f. Cooperative binding — uptake of clusters 


The binding of metallic elements to proteins could develop in a further 
way; in their higher oxidation states the hydrolysis of ions at pH 7 is quite 
extensive and some low molecular mass species of hydrolised species 
will be present in biological environments, e. g. surface waters. 


Inside cells (and also outside cells in the early history of life) sulphide 
radicals are also present and can react in exactly the same way as 
hydroxide radicals to give polysulphide units. 
We therefore predict and find polymeric clusters bound to proteins, such 
as Fe,O** (O* bridge) in hemerythrin, Fe,S, (S* bridges) and other Fe/S 
species in ferredoxins. 
The binding of groups of metal ions is opposed by entropic considerati- 
ons, but the unfavourable term in the association is greatly reduced if the 
anionic centres are already built into a pre-former, such as a protein, and 
it is further decreased if the protein folds in the adequate conformation for 
the uptake of the cluster. Different folds can favour different clusters, e.g. 
Fe,S, or Fe;S.. 

g. Transfer from aqueous media to proteins 


We have already comented on the possible changes on binding a metal 
or metal complex by a protein relatively to the same species in aqueous 
media. 

The results cannot be anticipated theorically, but the examination of 
Table VI shows the kind of effects that have been observed. 


Table VI 
Complexes or iron porphyrin IX with proteins 


protein binding E“/mV 
myoglobin (haemoglobin) His or His, HO >+100 
cytochrome b, His, His 0 
cytochrome Daas His, Met <+100 
peroxidase His or His, HO -200 
catalase Tyr or Tyr, HO -500 


As it can be seen, the range of redox potentials on binding the haem by pro- 
teins in different ways is more than 700 mV. Even without change of ligands, 
e.g. in the cases of peroxidase and myoglobin, the protein fold can distin- 
guish Fe(Il) from Fe(lll) by a differential factor of 10º. 

Similar results could be given for the series of iron/sulphur proteins in which 
differences in charge or even structure are amplified by the folding energy 
which, therefore, becomes also a selecting factor. 


An example helps to understand how the several factors operate to ensure 
the right distribution of required elements: that of blood plasma, where the 
metal ions are removed one by one by powerful ligands so as to avoid 
mutual interferences. Table VIl summarizes the main conditions and the 
basis of selectivity. 
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TABLE VII 
Carriers in blood plasma 


Metal ion Free M Carrier Basis of selectivity 

Iron (and manganese) 10º Transferrin Phenolates bind Mº 

Copper 10 An albumin lonised peptide >Nº binds Cu” 

Zinc 10º An albumin N'S donors in tetrahedral 
geometry bind Zn” 

Cobalt (as vitamin B,.) 10" Carrier protein The B, side chains are 
recognised 

Calcium 10º Phosphoproteins (when needed) Weak binding for most metais 

Magnesium, Sodium, 1040" None Freely circulating as lons 


Potassium,(Calcium) 

Thus, copper is removed first by albumin containing one or two ionised pep- 
tide bonds (>N') after complexation; although abundance is low, the geo- 
metry of the site (tetragonal) and the strength of the bond strongly favours 
copper relatively to zinc (the selectivity ratio can be as high as 10º). Iron is 
oxidised and the blood-protein transferrin picks-up Fe(lIll) using phenolate 
binding groups. Zinc can then be complexed by an albumin with three or 
4 nitrogen-donors in a tetrahedral geometry. Manganese may be transported 
by transferrin (same site as Fe(lll) but weaker bonds) and cobalt by cobal- 
min (vitamin B12). Sodium, potassium, magnesium and, to some extent, cal- 
cium are left free in the blood solution, but if extra calcium is necessary the 
solubility products of its carbonate or phosphate may be exceeded, hence 
there also is a special calcium carrier protein with phosphorylated side 
chains, generated under hormone control, that can transport this metal ion 
while maintaining the free Ca”* concentration at a sufficient low level, If all 
transition metals had not been removed first, these would bind strongly to 
the calcium sites and the process would be impossible. 


4. KINETIC EFFECTS 


Although exceedingly important, equilibrium binding effects are not the uni- 
que factors to be considered as determinants of the uptake of metal ions 
and their selection by biological systems. 


As mentioned before, kinetic considerations may restrict or favour the distri- 
bution of different elements into different cells and cell compartments. In 
general terms the following rates of various processes must be considered: 


a. In-flow of M to the environment of a cell, that may be assisted by certain 
ligands, e.g. siderophores for Fe” 


b. Out-flow of M from the environment of a cell, that may be assisted by 
other ligands, e.g. metallothionein to remove Cd** 


c. Rates of transfer in and out of cells involving formation of metal comple- 
xes and diffusion through membranes, pores and gates, which can be 
energy coupled,. 


d. Rates of (irreversible) binding of M by final sites, which can be assisted 


by insertion catalysts, e.g. the insertion of Mg”* in the chlorin ring of chlo- 
rophyill, 


The actual capture of M by a cell is no longer a thermodynamic problem as 
is frequently assumed to be. Rather we see it inserted in a steady state pro- 
cess in which uptake and rejection of wanted and unwanted elements must 
be considered simultaneously. 


The same can be said for the non-metals; their trapping must depend on 
membrane restrictions, thermodynamic effects and kinetics of several energy 
requiring reactions. Substances such as O,, N,. CO, or the ions SO |, NO; 
SeO ,, MoO, need to be activated and usually reduced before they can be 
bound into small covalent molecule intermediates, which may subsequently 
react or polymerise. Compounds such as Si(OH), and B(OH), are readly 
trapped by binding, probably to specific polysacharides. The remainning ani- 
onic species, e.g. F', Cl, and HPO ,, that cannot be oxidised or reduced are 
retained by membranes and controled by energised pumping in much the 
same way as Na* or Ca”. To phosphorus there is also left the possibility of 
condensation, forming one or two inevitably energived covalent bonds, e.g. 
R-OPO |, ATP, etc. It is clear, therefore, that the non-metal biological traps 
are based on kinetic stability to oxidation and hydrolysis, not on thermodyna- 
mic stability. 

Table VIll summarizes the situation concerning the use of thermodynamical 
and kinetical traps to capture and retain the various biological elements. 


TABLE VIII 
Traps for different biological elements 
Element Trap n 
Na*, K*, Mg”, Ca”, CI Physical retention of ions by membranes (fast chemical exchange) 


(H.PO,;, HPO,*, S0/) 
Cu*, Znº, Ni”, Mo, Si, Thermodynamically and kinetically stable binding to polymers (proteins 
Fe”, Mn” and polysacharides) 


Fe”, Co”, Ni, Mg” Thermodynamic and kinetic traps in small chelates, e.g. polypyrroles, 
bound to proteins 

HCNOP,S,Se,,Br Kinetic retention in thermodinamically unstable small or large molecules 
by covalent bonds 

SiCa,PC,S,F Precipitates (biominerals) 


5. CONCLUSION 


The above description of the selection of the chemical elements by living 
systems is based on present-day abundances and availabilities from the 
external environment. These availabilities, however, have changed dramati- 
cally from the earliest times to the present day partly due to the selections 
that life itself has made. It is appropriate to recall here that the evolution of 
an organism able to use light to decompose water generated dioxygen and 
this substance not only changed the availability of carbon, nitrogen and sul- 
phur but also altered that of the many metals, (e.g. iron, zinc and copper) 
and other non-metals. 


Life has, therefore, changed the atmosphere and the surface of the Earth, 
hence the composition of rivers, lakes, oceans and the local distribution of 
many elements in soils. Present availabilities are a reflexion of life's evolving 
steady state. 

lt is, then, to be expected that some types of early evolution will be found 
today only in residual ecological niches and these will nicely illustrate the 
theme of selection related to function. 


A good example would appear to be the association of nickel with 
CH,/CO,/H, chemistry in the archae-bacteria. Vanadium too has peculiar 
functions in a few ancient species: tunicates, toadstools (Amanita), brown 
and red algar, a lichen and a fan-worm. 


Evolution led to the replacement of these metals by more effective “modern” 
combinations, but was this the result of chance or of necessity due to pro- 
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gressive difficulties in obtaining the required elements whose availability was 
overtaken by others? The recent discovery of a vanadium nitrogenase gene 
that becomes expressed in conditions of molybdenum deficiency is sympto- 
matic, but we must also ask why nitrogen fixation is restricted to certain bac- 
teria, 


Evolution in general has led to the specialisation of most elements, particu- 
larly the metals, to optimise their functional value. A particular opportunity for 
new functions arose with organization and specialisation of cell compart- 
ments and cell types and proceeded as far as to separation among diferent 
species. Boron, for instance, which is essential for certain plants, appears to 
be totally absent in animals and the same happens, of course, with complex 
species such as magnesium chlorophyll. 


Evolution is then a many arrowed process driven by optimum functionality as 
determined by the diversity of the environment and by the evolved diversity 
of needs. Expressed in a different way, Life chemistry can not be regarded 
as a permanent set of reactions tending to a stable state of equilibrium with 
the environment, but as a dynamic process of continuous interaction of living 
species with the environment whereby the environment itself is changed and 
the changes are reflected back into the chemistry of the living species. 
Where we have living organisms there is a non-equilibrium situation at work, 
and therefore a kinetic approach to the selection of the elements is often 
more meaningful than a purely thermodynamic one. This point is frequently 
missed in studies e.g. of environmental polution, in which great efforts are 
put into monitoring concentrations of toxic substances without appreciating 
that the rates at which these penetrate cells or are rejected by cells are of 
overwhelming or indeed predominant importance. 
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RESUMO 


O presente artigo tem por 
objectivo, poder contribuir para 
um melhor conhecimento das 
estratégias seguidas pelas uni- 
dades de alojamento turistico, 
permitindo às diversas entida- 
des que podem ajudar a 
melhorar o serviço de quali- 
dade no alojamento turistico 
intervirem de uma forma pró- 
activa. 

Procedeu-se ao estudo da rela- 
ção dos seis grupos estratégi- 
cos de qualidade na Costa de 
Lisboa com os seguintes con- 
juntos de variaveis: (1) “desem- 
penho”, (2) “entidades que 
podem ajudar a melhorar o ser- 
viço de qualidade” e (3) “estra- 
tágias e politicas governamen- 
tais no sector do Turismo. 


COMPORTAMENTO DOS 
GRUPOS ESTRATÉGICOS 
DE QUALIDADE FACE A 
VARIÁVEIS DETERMINANTES 
NA INDÚSTRIA DE 
ALOJAMENTO TURÍSTICO 
NA COSTA DE LISBOA 


Maria Manuela Sarmento Coelho 


Instituto Supenor Técnico - Secção Autónoma de Economia e Gestão 
Original submetido em: 13 de Maio de 1996 


1. INTRODUÇÃO 


Os gestores das unidades de alojamento turístico actuam de forma diferente 
face às seguintes variáveis, consideradas importantes para o sucesso da 
unidade: 


1. Desempenho: 


1.1 Taxa anual média de ocupação de quartos; 
1.2 Taxa anual média de ocupação de camas. 


2. Entidades que podem ajudar a melhorar o serviço de qualidade: 


2.1 Governo central; 

2.2 Autarquias e Juntas de turismo; 

2.3 Associações de hóteis e similares; 

2.4 Associações de defesa do consumidor; 
2.5 Outras associações. 


3. Estratégias e políticas governamentais no sector do turismo: 


3.1 Importância atribuída pelos empresários; 
3.2 Consentâneas com a qualidade. 


No artigo sobre “Estratégia dos Grupos de Qualidade na Indústria do Aloja- 
mento Turistico na Costa de Lisboa” (1996) determinaram-se as estratégias 
implementadas pelas unidades de alojamento e constatou-se que há seis 
grupos. Estes grupos associam unidades que na sua gestão valorizam 
estratégias semelhantes. 


No presente trabalho pretende-se analisar a relação entre os grupos 
estratégicos de qualidade e as variáveis referidas. 


Para a realização do objectivo foi necessário realizar um estudo explorató- 
rio, que se concretizou através de um inquérito. 
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2. METODOLOGIA DE INVESTIGAÇÃO 


A população alvo foram as 357 unidades de alojamento turístico da Costa 
de Lisboa, que estão distribuídas em 26 subsegmentos, que correspondem 
a 12 segmentos da oferta conforme o quadro n.º 1: 


Tipo de | Tipo de alojamento Categoria Tipo de alojamento | Categoria | 
alojamento | 


Hotel: 5 estrelas Apartamento-Hotel 4 estrelas Estalagem: 5 estrelas 
4 estrelas 3 estrelas 4 estrelas 
3 estrelas 2 estrelas 
2 estrelas | Albergaria: 4 estrelas 
1 estrela Aldeamento Turístico: luxo 
1 classe Turismo Habitação 
3 estrelas 2“ classe 
2 estrelas Turismo Rural 
Apartamento Turístico: 1 classe 
Pensão: 4 estrelas 2.º classe Agroturismo 
3 estrelas 
2 estrelas Pousada 


Quadro n.º 1 — DOZE TIPOS DE OFERTA DE ALOJAMENTO TURÍSTICO 
Fonte: Alojamento Turístico — Guia Oficial da D.G.T. 


A dimensão da amostra para a Costa de Lisboa é de 54 unidades de aloja- 
mento, para o grau de confiança de 1=95% (este calculo foi apresentado no 
referido artigo na introdução). 


Obtivemos 63 respostas válidas o que corresponde a um grau de confiança 
de 1=96,76%. 

No caso vertente, o objectivo central da investigação, é a relação entre os 
grupos de estratégicos de qualidade, na indústria do alojamento turístico na 
Costa de Lisboa, e as variáveis determinantes. 

Para a análise quantitativa dos dados utilizou-se o programa “Statistical 
Programme for Social Sciences” — SPSS versão 6.0, recorrendo-se à téc- 
nica de análise de variância para determinar se os seis grupos são significa- 
tivamente diferentes no que respeita às variáveis “desempenho”, “entidades 
que podem ajudar a melhorar o serviço de qualidade” e “importância atri- 
buída as estratégias e políticas governamentais”. 

Neste sentido, testou-se a igualdade de médias de cada uma das variá- 
veis estratégicas nos 6 grupos, tendo por base a verificação dos seguin- 
tes pressupostos da análise de variância: 


1. a população tem uma distribuição normal, 
2. os grupos têm igual variância. 


Para a análise da homogeneidade de variância entre os grupos, testou-se 
a hipótese nula contra a hipotese alternativa, isto é: 


H,: As variâncias dos grupos são iguais. 
H,: As variâncias dos grupos são diferentes. 


e aplicou-se o método de Levene, considerando um grau de confiança 
A=95%, OU seja, para um nível de significância de 0=5%. 

Sempre que o nível de significância de Levene fôr superior ao nível de signi- 
ficância considerado, não há razões para rejeitar a hipótese nula, ou seja, 
há igualdade de variâncias entre os grupos. 

Seguidamente, realizou-se a análise de igualdade de médias entre os gru- 
pos, considerando a hipótese nula (as médias do grau de importância de 


cada factor são iguais nos seis grupos) contra a hipótese alternativa (as 
médias do grau de importância de cada factor são diferentes nos seis gru- 
pos), ou seja: 


Ho :M=H=H=M,=4=H, 
HMM UAU AUÇEM, 


Rejeita-se a hipótese nula, sempre que a probabilidade do teste F seja 
inferior ao nível de signifcância de 5%, o que implica que as médias dos 
grupos são diferentes. Também se efectuou o teste de comparação múltipla 
de médias de Scheffé por ser mais rigoroso que o teste F. O teste de 
Scheffé permite identificar os grupos cujas médias sejam estatisticamente 
diferentes. Os resultados da aplicação da análise de variancia simples One- 
way são apresentados no quadro n.º 2. 


Variável | teve Festa Eron, Teste de Scheflé 
ne F F 


DESEMPENHO 


Taxa média anual de ocupação 0,347 0,5678 | 0,7242 | Não há grupos significativamente 
“de quartos diferentes. 
Taxa média anual de ocupação 0,251 | 7369 0,1409 | Não há grupos significativamente 

| de camas na diferentes. 


ENTIDADES QUE PODEM AJUDAR A MELHORAR O SERVIÇO DE QUALIDADE 


Governo Central 0,277 49524 | 0,0008 | O grupo 3 é significativamente dife- 
rente dos grupos 4 e 5 | 

“Autarquias e Juntas de Turismo 0,066 3, 3886 0, 0098 | O grupo 3 é significativamente dife- 
rente dos grupos 2, 5 e 6 | 


Associações de Hoteis 0,724 1, EB09 0, 1553 Não ha grupos signilicativamente 
e Similares diferentes 


Associações de Defesa do 0 450 6 821 0 0001 O grupo 3 é diferente dos grupos grupos 
Consumidor 506 
| Outras Associações 0,648 2 8628 02500 | Não ha grupos significativamente 
| diferentes 


ESTRATÉGIAS E POLÍTICAS GOVERNAMENTAIS NO SECTOR DO TURISMO 


importância atribuída às 0,088 45581 | 0,0015 | O grupo 3 difere significativamente | 
“políticas governamentais | o do grupo 2 | 
Estratégia do governo 0,338 4 8166 | 0,1247 | Não há grupos significativamente 
consentânea com a qualidade diferentes 


Quadro n.º 2 - RESULTADOS DA ANÁLISE DE VARIÂNCIA ONEWAY PARA AS 
VARIÁVEIS DE “DESEMPENHO”, “DE ENTIDADES QUE PODEM AJUDAR A MELHO- 
RAR O SERVIÇO DE QUALIDADE" E DE “ESTRATÉGIAS E POLÍTICAS GOVERNAMEN- 
TAIS NO SECTOR DO TURISMO". 

Fonte: Inquérito de Qualidade à Indústria de Alojamento Turístico na Costa de Lisboa 


No quadro n.º 3 constam os valores médios das variáveis “desempenho das 
unidades”, “entidades que podem ajudar a melhorar o serviço de qualidade” 
e “estratégias e políticas governamentais no sector do turismo” relativa- 
mente aos seis grupos. 


pres 


20 


[te [gn [es | ços ónget tos [o 


Desempenho 


Taxa anual média de ocupação | 

de quartos 0,5278 0,4106 0,5376 0,4192 0,5106 | 0.4848 
Taxa anual média de ocupação 

de camas 0,5089 0,3711 0 5536 0,4188 0,5500 | 0,5405 


Entidades que podem ajudar a melhorar | serviço de qualidade (escala de 1 a 5) 


Ê Esociações Ra do 


* consumidor 3,8261 4,5000 20000 | 2.000 


| Importância atribuída pelos | 
empresários 3,7308 4,8889 3,1667 28000 | 3, 4206 | 


Consentâneas com a qualidade 30800 | 22000 | 26000 | 20000 | 20000 


Quadro n.º 3 - MÉDIAS DOS GRUPOS ESTRATÉGICOS POR VARIÁVEL DE “DESEM- 
PENHO”, DE “ENTIDADES QUE PODEM AJUDAR A MELHORAR A QUALIDADE”, E DE 
“ESTRATÉGIAS E POLÍTICAS GOVERNAMENTAIS NO SECTOR DO TURISMO" 

Fonte: Inquérito de Qualidade à Industria de Alojamento Turístico na Costa de Lisboa. 


4. CONCLUSÕES 


No artigo referido no ponto 1 concluiu-se que para um grau de confiança de 
1=96,76% é possível agrupar as empresas da Indústria de Alojamento Turístico 
da Costa de Lisboa com afinidades estratégicas de qualidade, em seis grupos: 


* Grupo 1 — Qualidade de tipo médio alto; 
* Grupo 2 — Qualidade de tipo baixo; 

* Grupo 3 — Qualidade de tipo bom; 

* Grupo 4 — Qualidade de tipo excelente; 

* Grupo 5 — Qualidade de tipo médio; 

* Grupo 6 — Qualidade de tipo médio baixo. 


Estes grupos foram confirmados pela análise discriminante e diferenciam-se, 
nomeadamente, quanto: 


* ao processo como efectuam as suas estratégias de gestão; 

* ao modo com se posicionam face à concorrência, no âmbito alargado da 
gestão global; 

* à satisfação dos clientes; 

* à pesquisa de mercado; 

* ao modo como efectuam a publicidade e a promoção dos seus produtos 
e serviços; 

* à forma como gerem a imagem de marca da unidade; 

* aos meios que utilizam para controlar e garantir a qualidade dos seus 
serviços e produtos; 

* ao planeamento da qualidade; 

* à avaliação dos resultados; 

* à dependência das centrais de reserva e dos operadores turisticos; 

* à gestão de recursos humanos, com ênfase na responsabilidade, na 
competência, no conhecimento, na participação, na formação e na satis- 
fação dos recursos humanos. 


No quadro n.º 4 apresentam-se as percentagens relativas ao total de unida- 
des do grupo e ao total da amostra. 
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Seguidamente focamos apenas as estratégias de qualidade mais relevantes 
de cada grupo, que foram quantificadas no referido artigo. 


* Grupo 1 —- Qualidade de tipo médio alto 
É o maior grupo e está distribuido por nove segmentos ou tipos de oferta 
de alojamento. 41,27% das unidades desta região turística pertencem a 
este grupo (quadro n.º 4). 


As principais estratégias utilizadas são: 

Pesquisa de mercado, comunicação da qualidade, atenção a gestão dos 
recursos humanos, gestão por objectivos, garantia de qualidade dos pro- 
dutos e serviços, satisfação dos clientes e posicionamento em relação à 
concorrência. 


* Grupo 2 - Qualidade de tipo baixo 
Está distribuido por três segmentos de oferta e representa 9,52% das 
unidades da Costa de Lisboa (quadro n.º 4). 
As suas estratégias mais relevantes são: 


Responsabilização, participação dos recursos humanos na gestão e 
dependência das centrais de reserva e dos operadores turísticos. 


* Grupo 3 —- Qualidade de tipo bom 
Representa 15,87% das unidades da Costa de Lisboa, correspondentes 
a cinco segmentos de oferta de alojamento. As pousadas e os hotéis 
apartamento de 4* pertencem na totalidade a este grupo (quadro n.º 4). 
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As suas estratégias principais são: 

Gestão por objectivos, incorporação da qualidade na estratégia global, 
satisfação dos clientes, posicionamento face à concorrência, controlo e 
garantia da qualidade, participação dos trabalhadores na gestão e sua res- 
ponsabilização, formação e treino dos recursos humanos, dependência das 
centrais de reserva e dos operadores turisticos e competitividade no preço 
versus qualidade dos produtos e serviços. 


Grupo 4 - Qualidade de tipo excelente 

Representa 12,7% das unidades da Costa de Lisboa. Este grupo elucida 
claramente que qualidade não é sinônimo de luxo, pois temos quatro seg- 
mentos de oferta: um hotel de 5*, três hotéis de 4*, um hotel de 2*, uma 
pensão de 4*, uma estalagem de 4* e uma albergaria (quadro n.º 4). 


A orientação estratégica é a seguinte: 


Gestão por objectivos, incorporação da qualidade na estratégia global, satis- 
fação dos clientes, posicionamento face à concorrência, controlo e garantia 
da qualidade, pesquisa de mercado e sua influência no plano de qualidade, 
comunicação da qualidade, qualidade como objectivo da gestão, publicidade 
da imagem de marca da unidade ou da cadeia, fomento da cultura da 
empresa, dependência das centrais de reserva e dos operadores turísticos. 


Grupo 5 — Qualidade de tipo médio 

E composto de três segmentos de oferta. Representa 9,52% das unidades 
da Costa de Lisboa (quadro n.º 4). 

As suas estratégias mais relevantes são: 


Qualidade como objectivo da gestão, pesquisa de mercado e sua influência 
no plano de qualidade, comunicação da qualidade, competitividade no 
preco versus qualidade do produto e serviço, 


Grupo 6 - Qualidade de tipo médio baixo 

Este grupo tem uma quota de 11,11% desta região turística, sendo consti- 
tuido por pensões (quadro n.º 4). 

As estratégias mais importantes são: 

Qualidade como objectivo da gestão, participação dos recursos humanos 
na gestão e sua responsabilização, formação e treino dos empregados, 
dependência das centrais de reserva e operadores turísticos. 


Seguidamente apresentamos as conclusões da análise de variância que permi- 
tiram determinar os grupos que são estatisticamente diferentes relativamente 
as variáveis “desempenho”, “entidades que podem ajudar a melhorar a quali- 
dade” e “estratégias e políticas governamentais no sector do turismo”. 


1. 


Desempenho 


Constata-se que, embora os grupos não sejam estatisticamente diferentes, 
há todavia diferenças entre eles, quanto às variáveis: 


1.1 Taxa anual média de ocupação de quartos: 
* O grupo 2 apresenta o valor mais baixo, ou seja 41,06%; 
* O grupo 3 apresenta o valor mais alto, ou seja 53,76%. 


1.2 Taxa anual média de ocupação de camas: 
“ O grupo 2 apresenta o valor mais baixo, ou seja 37,11%; 
o grupo 3 apresenta o valor mais alto, ou seja 55,36%, muito próximo 


do grupo 5 com 55,00%. 


Os valores médios são os seguintes: 


* taxa anual média de ocupação de quartos no triénio de 91/93 — 
48,18%; 


* taxa anual média de ocupação de camas no triénio de 91/93 — 
49,05%; 


O grupo 3 -— qualidade de tipo bom — é aquele que apresenta diferen- 
ças positivas maiores, em relação aos valores médios dos grupos, para 
as variáveis de desempenho consideradas. 


Este grupo tem cinco tipos de alojamento que oferecem uma qualidade 
de serviço boa. Para além disso, tem a segunda maior percentagem de 
integração em cadeias hoteleiras, o que lhe proporciona vantagens com- 
petitivas, nomeadamente políticas de comercialização mais agressivas, 
gestão disciplinada com controlo uniformizado da actividade e redução 
de custos globais. 


O grupo 4 - qualidade de tipo excelente — adopta a estratégia de um 
serviço personalizado e diferenciado, pelo que, em média, as unidades 
deste grupo necessitam de maior número de efectivos por quarto e, por- 
tanto, têm custos mais elevados com o pessoal. Poderia inferir-se que há 
menor desempenho, contudo, um serviço personalizado implica que o 
pessoal recrutado tenha maior variedade de especialização. 


O grupo 2 -— qualidade de tipo baixo — apresenta o valor mais baixo de 
custo anual médio com o pessoal, o que está de acordo com a sua 
estrutura de tipo familiar. 


- Entidades que podem ajudar a melhorar a qualidade 


Os seis grupos apresentam comportamentos diferentes, em relação às 
variáveis que podem ajudar a melhorar o serviço de qualidade, havendo 
diferenças entre eles. Algumas variáveis apresentam, no entanto, dife- 
renças significativas: 


2.1 Governo central: 
Apresenta diferenças significativas entre o grupo 3 e os grupos 4 e 5: 


* O grupo 5 tem o valor mais baixo (2,0), próximo do grupo 6 
(2,13). Para estes grupos, o governo central pode ajudar a melho- 
rar o serviço de qualidade num grau baixo. 

* Oo grupo 3 apresenta o valor mais alto (4,0) significando que para 
este grupo, o governo central pode ajudar a melhorar o serviço 
de qualidade, num grau alto a tender para elevado. 


2.2 Autarquias e juntas de turismo: 


Esta variável apresenta diferenças significativas entre o grupo 3 e os 
grupos 2,5 e 6. 


. 


o grupo 6 tem o valor mais baixo (3,14), muito próximo dos valo- 
res do grupo 2 (3,17) e do grupo 5 (3,20), concluindo-se que, 
para estes grupos, as autarquias e juntas de turismo podem aju- 
dar a melhorar o serviço de qualidade num grau médio. 

o grupo 3 apresenta o valor mais alto (4,8), significando que o 
poder local e as juntas de turismo podem ajudar a melhorar o 
serviço de qualidade, num grau alto a tender para elevado. 


2.3 Associações de hotéis e similares: 
Embora não haja grupos significativamente diferentes, constata-se que: 


* 


o grupo 6 tem o valor mais baixo (2,57), concluindo-se que, para 
este grupo as associações de hotéis pouco podem ajudar a 
melhorar o serviço de qualidade. 
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* O grupo 3 apresenta novamente o valor mais alto (4,2). Para este 
grupo o serviço de qualidade pode ser melhorado num grau alto, 
com a ajuda das associações de hotéis. 


2.4 Associações de defesa do consumidor: 


Também esta variável apresenta diferenças significativas entre o 

grupo 3 e os grupos 5 e 6: 

* os grupos 5 e 6 têm o valor mais baixo (2,00) pelo que, para 
estes grupos as associações de defesa do consumidor pouco 
podem ajudar a melhorar o serviço de qualidade. 

* O grupo 3 apresenta o valor mais alto (4,5), implicando que para 
este grupo estas associações podem ajudar a melhorar o serviço 
de qualidade, num grau alto a tender para elevado. 


2.5 Outras associações: 

Embora não haja grupos significativamente diferentes, constata-se 

que: 

* o grupo 6 tem o valor mais baixo (1,57), implicando que outras 
associações muito pouco podem ajudar a melhorar o serviço de 
qualidade. 

* O grupo 3 apresenta o valor mais alto (3,3). Para este grupo, 
apenas medianamente outras associações podem ajudar a 
melhorar o serviço de qualidade. 


Os valores médios são os seguintes: 


* Autarquias e Juntas de Turismo (3,90) — grau aproximadamente 
alto; 

* Associações de Hotéis e Similares (3,52) — grau médio a tender 
para alto; 

* — Associações Defesa do Consumidor (3,36) — grau médio; 

*— Governo (3,13) — grau médio; 

* Outras Associações (2,71) — grau baixo a tender para médio. 


O grupo 3 - qualidade de tipo bom — é mais uma vez aquele que 
sobressai, apresentando diferenças positivas maiores, relativamente aos 
valores médios dos grupos, para as variáveis “entidades que podem aju- 
dar a melhorar o serviço de qualidade”. 


Este grupo tem boas relações com o poder central e local e com “lob- 
bies”, o que lhos confere força politica e negocial. 


Neste contexto, a visão deste grupo baseia-se em que a unidade está 
inserida no meio empresarial, social e político e a intervenção das dife- 
rentes forças conduzirá a um melhor serviço de qualidade. Entendem 
que a qualidade não é apenas o serviço prestado pela unidade, mas 
todas as infra-estruturas que a rodeiam, que vão desde as estradas, aos 
parques de estacionamento, aos transportes, aos espectáculos, às expo- 
sições culturais, às diversões, aos restaurantes, aos bares, aos catés, às 
comunicações (telefónicas e correios), até aos hospitais e clínicas, pas- 
sando pela limpeza, segurança e estabilidade política. 

O grupo 1 — de tipo medio alto — posiciona-se a seguir ao grupo 3 e 
embora seja um grupo com 6 segmentos de oferta e representando 
41,27% da amostra, também perfilha a opinião de que a qualidade que 
se presta ao hóspede pode ser melhorada com a intervenção e com a 
colaboração das várias entidades interessadas no sector do Turismo, 
que actuam no meio envolvente às unidades. 


Note-se que o grupo 6 - qualidade de tipo médio baixo — apresenta 
neste conjunto de variaveis o valor mais baixo, sendo seguido pelo 


grupo 5 — qualidade de tipo médio. As unidades pertencentes a estes 
grupos praticam uma gestão essencialmente virada para o interior da 
unidade e pouco esperam dos grupos de influência e dos poderes cen- 
tral e local. 


Estratégias e políticas governamentais no sector do turismo 


As duas variáveis deste conjunto têm comportamentos diferentes em 
relação aos seis grupos. O grupo 3 é o que apresenta diferenças positi- 
vas maiores, em relação aos valores médios dos grupos, para as variá- 
veis de consideradas: 


3.1 Importância atribuída pelos empresários: 
O grupo 3 difere significativamente do grupo 2. 


* O grupo 2 tem o valor mais baixo (2,50). Os empresários deste 
grupo atribuem pouca importância às estratégias e políticas 
governamentais. 

* O grupo 3 apresenta o valor mais alto (4,89). Este grupo atribui 
as estratégias e políticas governamentais um grau de importância 
elevado. 


3.2 Consentâneas com a qualidade: 


Embora não haja grupos significativamente diferentes, constata-se 

que: 

* os grupos 4 e 5 têm o valor mais baixo (2,00), demonstrando 

que, para os empresários destes grupos, as estratégias e políti- 

cas governamentais são pouco consentâneas com a qualidade. 

* Oo grupo 1 apresenta o valor mais alto (3,08). Contudo, para este 
grupo, apenas medianamente as estratégias e políticas governa- 
mentais são consentâneas com a qualidade. 


Os valores médios são os seguintes: 


*— Importância atribuída pelos empresários (3,57) — grau médio a ten- 
der para alto; 

* | Consentânea com a qualidade no âmbito do Turismo e Alojamento 
Turístico (2,66) — grau baixo a tender para médio. 


Nas “estratégias e políticas governamentais no sector do Turismo” o 
grupo 3 - qualidade de tipo bom - e o grupo 1 — qualidade de tipo 
médio alto - são os que se evidenciam por terem os valores médios 
máximos. Para estes grupos a qualidade depende não só de condições 
endógenas, mas também das exógenas à unidade, condicionadas obvia- 
mente pelo poder local e central. 

O grupo 2 - qualidade de tipo baixo — pouca importância atribui às polí- 
ticas governamentais. 

O grupos 4 — qualidade de tipo excelente — e o grupo 5 — qualidade de 
tipo médio — apresentam valores médios baixos sobre “a estratégia do 
governo é consentânea com a qualidade”, na medida em que esperam 
mais das políticas governamentais, para poderem melhorar o seu negó- 
cio e consequentemente a economia nacional. 


A análise do sector permite-nos concluir que os empresários estão apreensi- 
vos quanto às macro políticas governamentais e especialmente os grupos 4 
e 5. Por ordem decrescente de importância são apresentados os motivos 
das preocupações quanto à: 


política sectorial: devido à necessidade de haver um ministro para um 
sector que é tão importante na economia nacional. De realçar que o 
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turismo contribui para o PIB em 5,9% (em 1994), emprega 250000 pes- 
soas (5% da população activa) e contribui em 50% das receitas cambiais 
na balança de transacções correntes. 

2. política fiscal: devido aos valores elevados da tributação do IVA (apesar 
de ter baixado para 12%), da Segurança Social e do IRC. 


3. política financeira: devido às taxas de juro elevadas e a alguma instabi- 
lidade quanto à inflação e política cambial; 

4. política de transportes e comunicações: porque as vias de comunica- 
ção, a sinalização e os transportes estão ainda aquém dos que outros 
destinos turísticos concorrênciais com o nosso têm; 

5. política de ambiente: porque é necessário melhorar a limpeza, as facili- 
dades e a segurança dos parques florestais, das reservas, da praias, dos 
rios e dos lagos. 


6. política de marketing: apesar do ICEP ter investido na publicitação de 


Portugal no estrangeiro é, porém, necessário continuar a trabalhar nesse 
sentido e dar uma imagem atractiva. 
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RESUMO 


Este trabalho descreve a pes- 
quisa objectiva orientada para 
a determinação dos processos 
de ajuste e centragem de com- 
ponentes ópticos. Os conceitos 
teóricos de ajustagem foram 
ordenados sequencialmente 
num sistema computacional 
CAA (Computer-Aided-Adjus- 
tment), para melhoria da quali- 
dade e obtenção de uma alta 
precisão mecânica no alinha- 
mento de lentes. 


ABSTRACT 


This work describes the objec- 
tive research oriented to find 
adjustment processes for cen- 
tring optical components. The 
theoretical adjustment concepts 
have been sequentially ordered 
in a computer system CAA 
(Computer-Aided-Adjustment), 
to improve quality and obtain 
high mechanical precision when 
straightening lens. 
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1. INTRODUÇÃO 


Os produtos de alta precisão da engenharia mecânica devem trabalhar com 
a mais elevada precisão. Para atingir estes propósitos podem ser consegui- 
das melhorias na qualidade pela utilização de medidas objectivas de ajusta- 
gem, quando se aplicam as conclusões obtidas através de um conceito de 
CAA (Computer-Aided-Adjustment) consistente, nomeadamente aplicado ao 
ajuste e centragem de mecanismos em aparelhos ópticos. 


2. OBJECTIVO DA CENTRAGEM 


Neste trabalho de investigação é tratado o problema da centragem automática 
de lentes em sistemas assistidos por programas computacionais. Para isto 
foram feitas pesquisas prévias sobre a teoria do ajustamento mecáânico e foi 
projectado o equipamento necessário. Além do estudo dos parâmetros que 
influenciam o processo, foram utilizadas as possibilidades de reconhecimento 
da centragem com o auxílio de Sensores de distância de Fibras Ópticas (SFO). 


3. CARACTERIZAÇÃO DA DESCENTRAGEM DE LENTES ESFÉRICAS 


Segundo a norma DIN 3140 existe a descentragem de uma lente quando o 
eixo de referência e a normal à superfície opticamente activa não se tocam 
no vértice de centragem, figura 1. 


O vértice de centragem é assim o ponto de uma superfície opticamente 
activa no qual ela é atravessada pelo eixo de referência. Uma lente esférica 
está desta forma descentrada quando o seu eixo óptico não coincide com o 
eixo de referência. O eixo óptico da lente é definido como a recta que passa 
pelos centros de curvatura das superfícies opticamente activas. O eixo de 
referência pode ser escolhido de várias maneiras diferentes. Se a superfície 
lateral da lente possui uma largura suficiente, o eixo deste cilindro lateral 
(eixo de forma) pode ser utilizado como eixo de referência; para cilindros 
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centro de curvatura 


eixo de referência 


eixo óptico 


vertice de centragem normal à superiidie 


Figura 1 Definição do erro de centragem 


laterais finos, ou para uma diferença muito grande na sensibilidade das 
superfícies da lente, é mais vantajoso definir-se como eixo de referência O 
eixo de simetria do contorno da superfície mais sensivel. Uma outra possibi- 
lidade é a selecção, como eixo de referência, do eixo de simetria do assento 
da lente (por exemplo uma secção anular). 


A medida para a descentragem é o ângulo «x entre a normal à superfície no 
vértice de centragem e o eixo de referência, figura 2. Este ângulo x tem 
valor igual ao ângulo de basculamento das superfícies activas em relação à 
posição ideal. Para raios de curvatura finitos das superfícies das lentes, x, 
pode ser definido, em graus, a partir do raio de curvatura r, e da distância 
perpendicular a, entre o centro de curvatura e o eixo de referência conforme 
a equação (1): 


Eixo de referência 
e eixo de rotação 


Figura 2 Medição do batimento da superficie por um Sensor de Fibras Ópticas (SFO). 
sin(x,)= 7. * E (i=1,2) (1) 
SO 60 


[1 
(F = factor de correcção de unidades) 


F =3438= 


AK = a, + 3438 (Nem Ab hem mm) (2) 


4. DETERMINAÇÃO DO ANGULO DE BASCULAMENTO DA 
SUPERFÍCIE DE UMA LENTE 


O processo para a determinação do ângulo de basculamento de uma super- 
ficie baseia-se no princípio da medição mecânica do batimento da superfície 
por meio de apalpação mecânica. Pela substituição do apalpador mecânico 


por um Sensor de distância de Fibras Ópticas (SFO), o processo de medi- 
ção é feito sem contacto e possui assim as vantagens de ser uma apalpa- 
ção sem atrito ou outras interferências mecânicas. 


A figura 2 mostra o arranjo para a medição do batimento da superfície com 
o auxílio de um SFO. A obtenção do ângulo máximo de basculamento x em 
graus é feita pela determinação do batimento Ab de uma lente que gira em 
torno do eixo de referência. 


A equação (2) vale para uma apalpação perpendicular à superfície da lente. 
O batimento medido da superfície, no caso de uma apalpação não-ortogo- 
nal, é sempre maior que o de uma ortogonal, isto é, o valor obtido aqui é o 
valor mínimo. 


5. SENSORES DE FIBRAS OPTICAS 


O princípio para a determinação da distância por meio de SFO consiste em 
captar e avaliar o fluxo luminoso que é emitido por um condutor de luz 
(LWL), reflectido por uma superfície e captado por um receptor. O funda- 
mento para isto é composto pela lei da medição fotométrica da distância e 
pela geometria do feixe luminoso. Quando dois núcleos, constituídos por um 
emissor e por um receptor de luz, são dispostos de tal maneira que as suas 
superfícies estejam na mesma direcção num plano ou paralelas a este, 
então considera-se ter um sensor de fibras ópticas (SFO) simples. Quando 
uma superfície reflectora penetra no cone de luz do LWL emissor, uma 
parte do feixe emitido é então reflectido. Se a distância entre esta superfície 
e o SFO aumenta, as superfícies dos cones de emissão e recepção come- 
çam a sobrepor-se a uma certa distância perpendicular (s.) do SFO, isto é, 
tem início a iluminação do LWL receptor (figura 3). 


Superfície receptora 
do LWL iluminada 


Cone de emissão 


supertície 
reflectora 


Figura 3 Superfície de um LWL receptor parcialmente iluminada 


Numa distância s, ocorre o pico óptico, quer dizer, é atingido o fluxo 
máximo de radiação & para uma iluminação total do LWL receptor. Na 
figura 4 está representada qualitativamente a dependência do fluxo de radia- 
ção coma, D=f(s). 
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S. - Inicio da sobreposição dos cones de emissão e recepção 


5, - Distância na qual o pico óptico & atingido 


Figura 4 Curva característica tipica de um sensor de fibras ópticas 


A faixa de interesse para as medições, com a grande inclinação, esta locali- 
zada entre os pontos s, e s,, especialmente onde a relação p=f(s) o 
quase linear. 

A resolução alcançável (na faixa de poucos um) e a faixa de medição são 
dependentes da construção do sensor óptico e das características de refle- 
xão da superficie. 


O grau de reflexão de superfícies de vidro, na faixa da luz visível, monta a 
cerca de 4%. Em superfícies não espelhadas este grau de reflexão pode ser 
inferior a 1%. Isto significa que para a obtenção de maiores resoluções 
devem ser utilizadas: uma intensa fonte de irradiação e um receptor sensi- 
vel. Por isto optou-se pela utilização de luz infravermelha como elemento de 
medição. 


6. REPRESENTAÇÃO ESQUEMÁTICA DO CONJUNTO 


Uma lente gira sobre uma casca anelar. Pela movimentação da lente, em 
relação ao eixo de rotação, sobre este anel, a sua superficie superior é cen- 
trada, A superficie inferior já está automaticamente centrada em relação ao 
eixo do anel, considerando-se que cada uma das superfícies da lente seja 
perfeitamente esférica. Um amplificador de emissão controla a pulsação de 
um LED. A radiação obtida é incorporada dentro de uma fibra óptica, dirigida 
ao ponto de medição e, por meio de um jogo de microlentes, focada sobre a 
superfície da lente. As fibras de recepção do SFO captam o feixe reflectido 
direccionam-o a um fototransistor. Um transformador-amplificador transforma 
o fluxo de luz, com a ajuda do fototransistor numa corrente, cuja amplitude 
representa uma medida para a distância entre o SFO e a superfície da lente 
e serve para a determinação da descentragem. Simultaneamente a posição 
angular do anel é determinada por um sistema de detecção angular, e, para 
cada incremento, uma amplitude de tensão é associada num sistema de 
aquisição de dados. O incremento ao qual foi associado o maior valor de 
tensão serve para a escolha de um dos dois motores deslocadores para a 
centragem da lente. Este motor é accionado e desloca a lente em relação 
ao eixo de rotação de tal forma que a distância entre o SFO e a superfície 
da lente neste ponto se torne menor. A descentragem da superfície da lente 
é desta forma reduzida. Este processo é repetido várias vezes. Caso a dife- 
rença de amplitude de cada incremento por rotação não se altere mais, ou 
caso um valor máximo da amplitude não possa mais ser determinado, então 


LO] 


a superfície da lente está centrada com a precisão necessária dentro das 
possibilidades do equipamento. 


SAD — Sistema de aquisição de dados 
RE — Receptor 

AE — Amplificador de emissão 

SOF — Sensor de fibra óptica 

CM — Controle do motor 

MD — Motor de deslocamento 

MAA — Motor de accionamento do anel 
ÃO — Anel de apoio 

CR — Controle da rotação 

SMA — Sistema de medição do ângulo 


Figura 5 Representação esquemática do conjunto 


7. ESTRUTURA CONSTRUTIVA DO EQUIPAMENTO DE AJUSTAGEM 


O equipamento de ajustagem está montado de tal forma que as direcções 
do movimento de ajustagem sejam perpendiculares entre si e na direcção 
do eixo de rotação. Como base para o movimento de ajustagem da lente foi 
prevista a rotação de um excêntrico. Para melhorar a movimentação da 
lente, foram colocados impulsores entre os excêntricos e a lente com um 
ângulo de inclinação que garante a fixação da lente sobre o anel. 
O esquema construtivo está mostrado no figura 6. 


Para uma transmissão sem folgas do movimento de ajustagem, estes 
impulsores são montados sobre molas de fita e podem ser pré-tensiona- 
dos pela regulação da mola de fita frente aos excêntricos. Uma outra pos- 
sibilidade de pré-tensão é dada pela regulação dos excêntricos (com para- 
fusos de ajuste). Estas possibilidades de ajustagem são necessárias para 
se poder adaptar o equipamento de ajustagem aos diferentes diâmetros 
das lentes. 
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Parafuso de ajuste 


Braço giratório 


Mola 
Impulsor 


Lente 


E 
CEEE TETE 


Exceêntrico 


Figura 6 Componentes de movimentação do equipamento de ajustagem 


A transmissão do movimento de ajustagem para a lente é feita por uma super- 
fície de contacto inclinada em cerca de 30º em relação ao eixo de rotação. 
Com isto consegue-se que e lente seja pressionada sobre o anel pela mola de 
fita, que age sobre o impulsor, com aproximadamente metade da força de fixa- 
ção. A disposição esquemática do excêntrico, impulsor e lente estã mostrada 
na figura 7. 


| Movimento de ajustagem 
| ses E 
4 FA Excêntrico 
ne | 
| 
A Avanço 
Z 
| FA Impulsor de ajustagem 
/ 
Anel Z Mola de fita 
de apoio / 
| iii Z Base deslocável 
| E mm 


Figura 7 Representação esquemática da unidade de ajustagem 


8. O PROCESSO DE AJUSTAGEM 


A ajustagem é possível tanto manual quanto automaticamente. Antes da ajus- 
tagem devem ser definidos os critérios de interrupção do processo necessá- 
rios, Estes critérios de interrupção são o valor mínimo permissível da tensão 
para a amplitude do batimento da superfície Aymin [Nm] e o número máximo de 


revoluções Nmaz - O valor da tensão é tomado pelo SFO em [mV] e transfor- 
mado pelo computador em [um]. Com Aymin determina-se a precisão da cen- 
tragem. 

Um algoritmo de ajustagem está mostrado no figura 8. A distância SFO- 
lente é registada em forma de uma amplitude de tensão Aymax pelo SFO 
durante uma rotação do anel de suporte. Se Ay, foi determinado muito 
pequeno, o programa não pode interromper o processamento, porque 
pequenos impulsos perturbadores que agem sobre a tensão de emissão tor- 
nam o valor da tensão determinado AYa« Sempre maiores que AYmin: À lente 
seria desta forma continuamente deslocada em torno do vértice de centra- 
gem. Por isto existe a possibilidade de se limitar o numero de revoluções 
admissível. O processo de ajustagem é então interrompido quando se atinge 
o número máximo de rotações. 


Inicio 


Excêntrico na posição inicial 


determinar Aya Na 


| Determinar Aya pOr rotação 
Centragem 
| Sm 


Seleccionar motor 
| Passo de deslocamento 


Contar rotações sim 


n “ mai 


| | não 


A partir da distância SFO-lente é produzida uma amplitude de tensão, pelo 
SFO, durante uma rotação do anel de suporte. Simultaneamente, a cada 
incremento da rotação é associado um valor da tensão. Estes pares de 
valores são armazenados. No fim de uma rotação os valores máximo e 
mínimo da tensão são verificados e obtém-se a diferença AY max 


Figura 8 Algoritmo de ajustagem 


A coordenação de um dos dois motores e a determinação do dispositivo de 
deslocamento é feita de acordo com a posição da maior amplitude e respec- 
tivo ângulo, o qual é obtido por meio de um disco graduado. Esta posição 
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marca a localização da maior distância entre o SFO e a lente. A determina- 
ção sobre qual dos motores será accionado é feita de acordo com as 
regiões fixadas na figura 9, 


90 o 0º 
Disco graduado 


(22,5) Fá 


Figura 9 Regiões para controle do dispositivo de ajustagem 


O equipamento permite obter 90 incrementos por rotação do disco gra- 
duado. A medição e a distribuição dos incrementos são feitas sempre a par- 
tir de um mesmo ponto inicial. Os limites das regiões, com respeito a este 
ponto de partida, são determinados como a seguir: 


Região Mol: Incremento 1 - 22 (0º - 90º) 
Região M,R: Incremento 23 - 45 (90º - 180º) 
Região M,R: Incremento 46 - 68 (180º - 270º) 
Região ML: Incremento 69 - 90 (270º - 360º) 


Por exemplo caso o valor máximo associado ao incremento nº 36, é accio- 
nado o motor 1. Vendo-se o processo a partir do motor 1, a lente sofre um 
pequeno deslocamento na direcção do eixo de rotação (direcção L). 


Após o término deste passo de ajustagem, n é comparado com Nmax-: Se n > 
Nmax: O processo é interrompido. Caso não tenha sido interrompido, a tole- 
rância da ajustagem AYmex é novamente detectada e verificada. Se AYmax É 
maior que AYmin, é feita novamente uma ajustagem. Caso contrário, o pro- 
cesso é interrompido, isto é, a lente está centrada. 
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ABSTRACT 


Existing methods of calculation 
fire resistance of timber mem- 
bers in walls are based on the 
sum of the contribution indices 
of its various components, 
(linings, claddings, sheathings 
and timber studs) to the stabi- 
lity, integrity and insulation per- 
formance of the whole element, 
Older forms of construction, 
such as wooden frame and 
stone masonry walls and timber 
lath and plaster walls are not 
commonly considered due to 
the variation of their condition 
which could render any evalua- 
tion worthless. 

However in some situations of 
city rehabilitation these ancient 
construction systems are pre- 
served and similar walls using 
new materials with traditional 
techniques are built in a spirit 
of structural conservation (Pie- 
tures 1 - 2 Appendix), 

It is therefore important to esta- 
blish some classification criteria 
to permit that old systems are 
to compete in equal terms as 
new ones in relation with fire 
regulations and thal some 
sense of security could be 
given to historic centre's inhabi- 
tants. 

The present repor! was based 
on an experimental study in 
which a masonry wooden frame 
wall section, plastered in both 
sides with mortar, was tested in 
an ISO 834 time dependent 
temperature evolution until its 
full destruction. Temperature 
measurements and relevant 
aspects were registered and 
results were correlated with a 
heat transfer analytical model. 
The obtained results showed 
that an effective temperature 
stop lining was achieve with the 
mortar layer and the reduction 
of the wood section to approxi- 
mately one half, inside the wall 
section, only occurred after 2h. 


FIRE RISK IN 


HISTORIC BUILDINGS 


José Gonçalves 
Eng. M in Lisbon Council 
Original submetido em: 24 de Junho de 1896 


Contents 


1. Theory of wood combustion and mortar behaviour in high temperatures; 
2. Analytical model used to simulate the tests; 

3. Experimental tests; 

4. Correlation of results; 

5. Conclusions; 


1. THEORY OF WOOD COMBUSTION AND MORTAR BEHAVIOUR IN 
HIGH TEMPERATURES 


Existing analytical methods to simulate the stress-strain behaviour of wood 
when exposed to fire are based on calculation of the time to ignition, rate of 
charring, heat and moisture transfer and mechanical behaviour in loadbea- 
ring condition. Ignition can occur in two different ways, either when the flam- 
mable mixture of decomposition products from wood and air, reaches the 
flame (pilot ignition) or when the flame is in direct contact with the surface 
(surface ignition). However the parameters studied in ignition modelling were 
basically the same, ie: 


— Ignition time; 

— Surface temperature; 

— Temperature distribution inside the material; 
— Loss of weight; 

— Absorption of radiant energy; 

— Moisture ratio, and 

— Convective effects. 


Models by Simms [1]. [2] showed that the intensity of radiation at which igni- 
tion is just possible, might expected to vary between 0.25 and 0.33 cal em * 
s being 0.30 cal cm * s ! the value normally considered for practical calcu- 
lations. For this intensity of radiation the ignition time before collapse should 
be at least 10 min for thick specimens, providing some time for fire control. 
Wesson et all [3], correlating the spectral absorptance relative to the incident 
radiation as a dominant variable for ignition of dry wood, concluded that the 
ignition time is approximately proportional to the wood density and a single 
correlation sufficed to predict the pilot ignition time of all wood species. The 
same conclusion was reported by Koohyar [4] who studied the spontaneous 
ignition of wood with specimens of pine, fir, mahogany and redwood. 


Models for wood pyrolysis showed the existence of four characteristic zones, 
parallel to the heated surface in which the basic physical and chemical 
mechanisms are expected to occur. These mechanisms can be summarise 
as internal convective heat flow, thermal conduction, storage of heat, 
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exchange of heat on material surface, vaporisation of water and chemical 
reactions inside the wood. The characteristic zones are illustrated in Figure 
1.1 and each one correspond to a temperature level at which the above fac- 
tors change in accordance to the biological variation of wood. 


N 
O 
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Jo 


Zone € 


Zone D 


Original 


Surface 


Figure 1.1 Characteristic zones for pyrolysis of wood [5]. 


Zone À is associated with temperatures below 200ºC where the gases for- 
med during pyrolysis process are not ignitable but a slow loss of weight and, 
in certain cases, a char surface begins to occur. 


Zone B corresponds to a temperature range between 200 and 280ºC where 
the flow of gases both towards the surface and further into the wood became 
noticeable. The temperature at which the net result of all pyrolytic reactions 
becomes measurable has been taken as one of the definitions of the ignition 
point of wood. 


in Zone C, to temperatures between 280 and 500ºC, a complex process of 
heat and mass transfer occur due to the large emission of COZ and water 
vapour. A set of cracks of a depth which increases in step with the expan- 
sion of the pyrolytic zone, develops at the surface and exchange of radiant 
and convective heat between the external heat source and the cracked zone 
can now be transferred directly to the uncharred wood by way of these 
cracks. At this stage combustion with a flame occurs in the gas phase, out- 
side the wood, due to the rate increase of emitted gases. The charcoal layer, 
due to its lower thermal conductivity, retards penetration of heat and the 
occurrence of an exothermic state further inside the wood. 


In Zone D, to temperatures above 500ºC, the charcoal glows and is oxidised 
until it reaches a temperature over 1100ºC. At this temperature its consum- 
mation is so fast as the penetration of the reactive zone. When all the wood 
is consumed, the luminous diffusion flame cease and only the non-luminous 
flame of the burning carbon monoxide and hydrogen remain. 


Models purposed by Kansa [6], Alves e Figueiredo [7] and Roberts [8] assu- 
mes a multi-step pyrolysis reaction scheme based on the fact that all physi- 
cal properties of the wood, i.e., heat capacity, density, thermal conductivity, 
porosity, permeability and molecular weight are linear functions of the local 
char or the fuel concentration and there is a thermal equilibrium between 
solid and volatiles within the wood. 


Temperature values for constituents decomposition are obtained and some con- 
clusions show that ignition of wood is controlled by the behaviour of the most 
reactive component (hemicellulose) while the time require to complete pyrolysis, 
at a given temperature, is controlled by the least reactive (lignin). Further work 
by Simms [9] with Oak, Colombian pine and European whitewood show that 


of testing, which is a good indi- 
cator of the fire resistance that 
can be expected in similar 
constructions. 

Temperature measurements on 
the unexposed lining and the 
rate of smoke liberated also 
showed that a fire performance 
of 30:180 min. (integrity and 
insulation resistance), could be 
anticipated if the experiment 
had followed the fire testing 
standards which was impossi- 
ble to achieve due to laboratory 
restrains. 


the range of intensities at which ignition occurs is strongly dependent of the 


moisture content. Values for pilot ignition are show in Table 1.1. 


(Quercus sp.) 


| Colombian pine 
(Pseudotsuga 


taxifolia 


European 
whitwood 
(Picea abies) 


Moisture 
content 
percent 


Range of 
intensities 
cal cms” 


0.38-0.50 
0.55 

0.59-0.65 

0.69-0.75 


0.38-0.40 
0,50-0,55 


“0,41-0,45 


0,46-0,50 
0,54-0,55 
0,63-0,70 


| 0,74-0,76 


0,39-0,40 
0,43-0,45 
0,39-0,41 
0,45-0,50 
0,52-0,60 
0,57-060 
0,67-0,70 


0,39-0,39 
0,40-0,52 


ignition 
times, s 
415-140 
605 
635-530 
510-435 
1260-1115 
1020-630 
2580-2020 
430-160 
460-500 
380-180 


| St0-140 


2130-1440 
1940-1770 
3540-2230 


240-180 
610-370 
550-260 
3550-270 
2380-1520 
1800-300 


0,42-0,50 1520-530 
Table 1.1 Results obtained by Simms for pilot ignition. 


From the present data is possible to see that, taking as an example the 
Colombian pine, and for the same density and thickness, the increase of the 
moisture content in 60% is related with an increase of the intensity values to 
more than 0.70 cal cms”, clearly higher than the mininum necessary to 
signition (0.30 cal cms”). This proof the importance of this factor to the igni- 
tion of wood. 

Concerning the structural properties of wood, papers by Fu-Kuo Chang [10] 
and [11], related stress and mechanical behaviour of wood with time exposi- 
tion to elevated temperatures, concluding that at temperatures of about 
150ºC most specimens collapsed within five minutes of exposure while for 
temperatures of about 300ºC all the loaded specimens failed within 60 
seconds of exposure. According to the work of Schaffer [12], this is due to 
the depolymerization of cellulose with duration of heating for temperatures 
above 288'C.Therefore this temperature has been taken as a limit for the 
structural use of timber after a fire. 


In relation with the behaviour of cement based mortars when subject to ther- 
mal and structural loading, Piasta [13] investigated the decomposition of 
cement pastes under thermal exposure in a temperature range of 20 to 
800ºC and suggested that the shrinkage process is associated with two pha- 
ses regime. The first, for temperatures up to 100-250ºC, corresponds to an 
expansion of strains especially in a fast rate of heating while the second, 
associated with temperatures of 250 up to 670ºC, strains show contraction 
mainly due to the dissociation of Ca(OH).. At temperatures close to 500ºC 
there is an increase on formation of microcracks while, for temperatures 
close to 800ºC, the total porosity raise to values up to 0.223 emº/g. This irre- 
versible characteristics suggest that the structural use of cement based mor- 
tars is limited to temperatures below 670-700ºC. 
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2. ANALYTICAL MODEL USED TO SIMULATE THE TESTS 


The analytical model was based on the Fluent computer program [14] which 
simulate the combustion process inside the furnace and the heat transfer 
process within different solid regions. This is achieve by solving the conser- 
vation equations of mass, momentum, energy and chemical species using a 
simple control volume based difference method. 


The fundamental equations for conservation of energy in a conducting solid, 
including heat flux and volumetric heat source can be describe as: 


dT d dT d dT 
O » — cm of oo —s «jo E. = 0 
DC de TR DR COR (1) 


Where: 


T - temperature, “K,; 

t - time, S; 

p - density, kg mº; 

c - specific heat capacity, J kg'K”; 

K - thermal conductivity, Wm'K'; 

q” - volumetric heat source, Wm*; 

x,y - coordinates in the x and y directions, respectively, m. 


The equation for conservation of energy in the form of a transport equation 
for the static enthalpy h in the fluid, is as follows: 


q q o dT d D du 
-— ph-—(puh)+—(k+k)—-— 3 h JJ +4+—L+47,—L+8, =0 2 
dt i x, (4 od, o, dx : PIAS D ! dx, h ( ) 


With: 
h = Emis E ai fe ,dT (3) 


Where: 


h';,- partial enthalpy of formation of species /' 

m - mass fraction of species /” 

C,, - specific heat at constant pressure of species /' 
T,.- reference temperature; 

u - velocity, m s”; 

k - molecular conductivity, W m'K''; 

k, =u,/Pr, - effective turbulent conductivity, W m'K”; 

J - diffusion flux of species j'; 

Tu - VISCOUS Stress tensor; 

Sy Source term, which includes the heat of chemical reaction, (S,, reaction) and 
the exchange of heat with the dispersed second phase: 


no 
Simm E + fe, a R 
(4) 
Where: 
R, volumetric rate of creation of species j', 


M,- molecular weight of species ). 
hº, - enthalpy of formation of species j” 


The reported enthalpy value, h' is then computed as: 
h =3m (h, +h)] (3) 


|) 


E | 
=5m, fe, dT+ho IM, + Jem dT (6) 
) T Tr 


Where T,.,, is the reference temperature at which hº, is defined. 


Turbulent analysis based on the (k-£) model and using the Boussinesq hypo- 
thesis to calculate the Reynolds stresses were used in order to obtain an 
ensemble-averaged momentum equation for the fluid: 


d = dp dr, o dp—-— 
(Pt) + Pt; )= Cr AM F ET o dé ,) (7) 


Where: 


T; - VÍSCOUS stress tensor; 

p - static pressure; 

q, - gravitational acceleration in the | direction; 

F, - external body forces which arising as a result of the interaction with the 
dispersed phase on the fluid in the / direction; 


(pu, U,) - “Reynolds stresses” representing the effect of turbulence in the 
time-averaged (or mean-flow) values. 


Radiation was modelling using the Discrete Transfer Radiation Model 
(DTRM) where the radiation intensity is calculate for each position, s, and 
solid angle, w, as: 


j 


a + fe 
É + + [ Hs. uJdw (8) 
K 4x L 


dio. 
——+(a+o )l=a: 
ds 


Where: 


| - radiation intensity from the flame; 

a - absorption coefficient; 

o. - scattering coefficient; 

o - Stefan-Boltzmann constant (5.672 x 10 *“W/mK'); 
(a+0.) - optical thickness or opacity of the flame; 


Integration of (8) along a path for all wavelengths was executed and incident 
radiation heat flux was calculated assuming that all surfaces were diffuse. 
Two layers near the wall region were assumed considering a laminar viscous 
sublayer near the wall and a turbulent log-law sublayer at some distance of 
the wall. Discussion of the theoretical model is develop elsewhere [14]. The 
geometric grid used in the model corresponded approximately to a 2x2 em 
cell in the experimental apparatus. Near the wall model and inside it, where 
differences in material occurred, a higher definition was used. This grid is 
used by the analytical model to calculate the relevant parameters as tempe- 
rature, velocity vectors and pressure associated to the cells of the model. 
Types of cells summarised in Table 2.1 define the way how the cell is trea- 
ted during the solution process. 
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Cell Type Description Designation | Application 


INLET Positive Velocity Boundary Condition E 
(VELOCITY)-1 

INLET Negative Velocity Boundary Condition Combustion 
products 
Outlet from 
the Domain 


(VELOCITY)-2 

CONDUCTING | Wall Boundary Condition — Type 1 (Heat Mortar layer 
WALL Transfer Solved via a Conducting Equation) 

CONDUCTING | Wall Boundary Condition — Type 2 (Heat 

WALL-2 Transfer Solved via a Conducting Equation) Wood frame 
CONDUCTING | Wall Boundary Condition — Type 3 (Heat 

WALL-3 Transfer Solved via a Conducting Equation) 

WALL Wall Boundary Condition Combustion 


Chamber 
LIVE (FLUID) | Fluid E 


Wall 
Table 2.1 - Cell types and labels used in the model. 


Wi 
W2 
W3 


Gas and 
Combustion 
Products 
inside the 
Combustion 
Chamber 


3. EXPERIMENTAL TESTS 


In order to verify the integrity and insulation characteristics of a ancient woo- 
den wall, two 80x80x16 cm real scale models were executed, like is show in 
Picture 3. 


The materials and the procedure used in model's preparation were the same 
as those used in reconstruction of ancient buildings, and are summarised as 
follows: 


A — Structural frame of pine wood with a mean density of 500 kg/m3. 
B - Stone assembled with cement-lime based mortar, composed by: 


1 volumetric part of Portland Cement; 
1 volumetric part of hydrated lime; 
5 volumetric parts of sand. 


C - 3 mortar layers applied in the following order: 


C1 — First layer with the same proportion of the mortar described 
above, 5 mm depth. 


C2 - Second layer, 15 mm depth, with the following composition: 
1 volumetric part of Portland Cement; 
1 volumetric part of hydrated lime; 
6 volumetric parts of river sand; 


C3 — Final layer, 10 mm depth, with the following composition: 


1 volumetric part of Portland Cement; 
2 volumetric parts of hydrated lime; 
9 volumetric parts of river sand; 


B 


Picture 3 Model construction 


In order to achieve good adhesion among the three layers (C1-C2-C3), some 
time between preparation of each one was followed, according to the proce- 
dure normally observed: 


— Between layers C2 and C1: 3 days; 
— Between layers (3 and C2: 6 days. 


Before execution of this layer, a chicken wire mesh was nailed to the woo- 
den frame to improve the cement adhesion and this was very important to 
prevent premature collapsing in a fire situation. 


Model's geometry is show in Figure 3.1. Inconel sheathed type K thermocou- 
ples positioning, inside the wall section, are also visible in this Figure. 

To achieve some information about the effect of the lime mortar, other model 
was made with a mortar lining composed with only cement and sand. 
Results showed that an interesting effect in retarding the temperature raising 
was achieved with the lime mortar model which make this material a better 
insulator. 


The test facility consisted in a 115x115x115 cm internal section furnace 
connected to a natural gas burner with two main extraction ducts and one 
refrigeration system in order to cool the metal walls during test. The 
model's positioning inside the furnace was crucial due to various factors 
like heat exposition, flame propagation, surface cracking and structural 
integrity. Standard procedures require also a model positioning during 
tests as close as possible to the reality due to side effects that can be 
expected. Additionally, in order to prevent the heat transmission throug- 
hout the sides of the model surface to the back side of the wall, all the 
model was insulated with ceramic fire bricks and fire wool. With this pro- 
cedure the temperature readings were in fact those from the wall section 
and not those from the heat transmited to the furnace walls. Thermocou- 
ples disposition in furnace apparatus and inside the wall section are show 
in Figure 3.2. 
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qd 


SECTION 1 


0.03 


0.80 


SECTION 3 


0.030.100.03 


0.16 


x - Thermocouples 


Figure 3.1 Geometry of the models and thermocouples positioning. 
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Figure 3.2 Thermocouples numbering and disposition inside the furnace apparalus. 


Significant occurrences during test evolution are show in Table 3.1. The 
main factor is related with the COZ production outside the combustion cham- 
ber after 2h. testing. This was the cause for the test interruption before signi- 
ficant temperatures had been achieve on the unexposed side of the wall. 


Time | Temperature Occurence 
(min) | | é é | 


a Ti First Rs of smoke outside a Tramber 


59 943 Signals of resin from the wood combustion inside = 
chamber o o | 
71 968 | Wood combustion continues. Signals of resin outside the 
| chamber, 
| External layer starts to collapse. 
B8 1010 Smoke production continues outside the chamber. More 
| resin outside the chamber. 
| 116 | 1059 | | Long crack (20 cm aprox.) appear on model's surface 
119 "1068 | Intense smoke outside the chamber. Glowing gases on fur- 
| | nace chimney 
138 | Excessive production of smoke from the furnace to the 
exterior. 
Model's surface ignition. Test stop. 


Table 3.1 Significant occurrences in Test 1. 


However in Test 2 the full temperature regime was obtained after 6h. of con- 
tinuos burning and the temperature evolution for thermocouples disposition 
describe above, is show in Figure 3.8. 


From temperature evolution in thermocouple N.8 is possible to see a good 
agreement with the normalise temperature evolution of ISO 834. 


Thermocouple N.12 represents the temperature between the mortar layer 
and the wooden frame and is visible that a significant insulation effect was 
achieved by this layer. In fact to temperatures of 200ºC the insulation time of 
the protective layer can be estimated as 20 minutes, but to temperatures of 
400ºC the insulation time of this layer raises to 35 minutes and to temperatu- 
res of 600ºC the insulation time obtain is 60 minutes. Moreover the tempera- 
ture only reach 700ºC near the wood region after 150 minutes of testing and 
the value of 800ºC is only achieve after 285 minutes which is a good indica- 
tor of the insulation confered by the exterior layer. 


If the temperature inside the wood section is consider, we can see a slow 
heating process to temperatures up to 140ºC for the first 3h. period and then 
a considerable raise to temperatures up to 727ºC, after 6h. This is explain 
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by the slow loss of weight associated with the flow of gases towards the 
exterior and further into the wood section as we have seen in Figure 1.1, to 
temperatures up to 200"C. The raise in temperatures after 3h. of testing is 
associated to the emission of CO2 and water vapour by the cracks which 
appear at temperatures close to 500ºC. 


Experimental test 


Figure 3.3 Temperature evolution in Test N.2 


If the protective layer cracks, as it's visible in Picture 4, then combustion with 
a flame occurs outside the wall due to the rate increase of emitted gases 
from the wood section. This was visible in Test 1 after 138 minutes (Table 
3.1). Temperature raising in thermocouple N. 13 after 3h. exposition also 
confirm that fact. 


An interesting temperature evolution was achieved with thermocouple N.14, 
which was positioned inside the rock section of the wall. In this case an 
almost linear temperature evolution was obtained as we can see in Figure 
3.3 and this is explain by the higher thermal conductivity of the rock when 
compared with wood to temperatures up to 200ºC. To higher temperatures 
however the results show a much more steady behaviour for the rock due to 
the stability of it's physical properties. 

Another interesting temperature evolution was achieve with thermocouple 
N.15, placed between the wooden section and the unexposed side of mortar 
layer. In this case an almost constant temperature regime during the first 5h. 
period was obtained due to the lower conductivity of the wooden section until 
some wood still exists inside the wall. After it's consummation is expected 
that the temperature will rise like what happen in the central section. This 
was not obtained due to the establish limit for test term being 6h. 


Nevertheless is worth to see that a higher temperature value in the back 
section of the wall is always obtain. If we consider that an effective insulation 
was achieved with the ceramic wool and the fire bricks placed between the 
model sides and the furnace walls then we can conclude that the wooden 
section acted, not only as a structural element for the first 3h. of fire, but 
also as an effective temperature barrier to the remaining section of the wall. 


Temperatures on the back side of the wall were measured in two different 
forms, with a single point (thermocouple N.9) and taking an average value in 
5 different positions. This was intended to establish some kind of credibility 
in measuring temperatures for classification purposes. 


Picture 4 Cracks on exposed side of the model after 138 min of testing 


In fact, if we suppose that the temperature reached 140ºC after 220 minutes, 
we can consider this wall as a thermal insulation element of 180 min. and, 
with the results of Table 3.1 for smoke production, we can assign the follo- 
wing fire performance to this type of element: 


(30:180) 
(integrity and insulation resistance's), 


The problem which arises from the present analysis is related with the stan- 
dard procedures that would have been necessary to have follow in order to 
attribute a fire resistance classification. This was not possible to overcome 
due to laboratory circumstances 


4. CORRELATION OF RESULTS 


The obtained results from experimental analysis were correlated with the 
theoretical model in order to validate the assumptions made during the test 
procedures. 


The values for materials properties, used in analytical model, were taken 
from references, namely [15], [5] and [16] in relation to the mortar, wood and 
rock properties, respectively. 
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For mortar it was assumed a temperature dependant thermal conductivity 
with values varvying from 0.597 to 1.2548 W m' K' for temperatures varying 
from 293 to 1487 K. The specific heat was assumed a constant with a value 
of 862 J kg! K! and density with a value of 1990 kg mº. 

For wood the thermal conductivity was assumed as temperature dependent 
with values varying from 0.1250 to 0.1630W m” K” to temperatures varying 
from 273 to 673 K. For temperatures over 1073K, a constant value of 
0.052W m” K! was assumed as the thermal conductivity of char. In relation 
with the heat generation it was assumed a time dependent model, according 
to Table 4.1. 


Heat generation 
(Wim) 


7.60E+03 


2.20E+03 
2.00E+04 
2 .50E+04 
2.50E+04 


Table 4.1 Heat generation of wood for temperatures of 284-542ºC. 


Finally for rock it was assumed a temperature dependant thermal conductivity 
with values varying from 0.934 to 1.4424 W m'! K' for temperatures varying 
from 293 to 1487 K. The specific heat was assumed a constant with a value 
of 837 J kg' K! and density also a constant with a value of 2308 kg mº. 


Obtained results from both analytical and experimental models are show in 
Figure 4.1. 


Temperature evolution in experimental (T) and analytical (tp simulation 


Figure 4.1 Correlation of results of experimental and analytical models 


From this graphic its possible to observe that a relatively good correlation is 
achieve with analytical modelling although for the back side of the model 
values obtained showed a poor agreement since the begining of simulation. 
This was mainly due to the limitations of the model to simulate temperature 
dependency of properties in solid materials. However this constraints, a rea- 
sonably correlation for temperature inside the wood section was obtained for 
temperatures below the pyrolysis region (thermocouple N.13) and is possible 
to validate the results of experiments, with the assumed properties definitions. 


Moreover an interesting factor related with the apparent inversion of ten- 
dency in wood regions for temperatures close to the ignition point seems to 
occur between analytical and experimental data (thermocouples N.s 13 and 
14). This is attributable to the nature of the analytical model which is too sta- 
tic for wood combustion simulation above this temperature regime. In stone 
region, oppositely, a fairly agreement was obtained specially for temperatu- 
res above 300ºC (thermocouple N. 15) and for mortar region a very good 
agreement was achieved for all temperatures (thermocouple N.12). 


5. CONCLUSIONS 


The results of the tests are a good indicator for the use of wooden structural 
walls in historic buildings by two main reasons. 


The first is related with the considerably long time before occurence of tem- 
peratures higher then 140-180ºC on the unexposed side of the wall. This is 
a proof of its considerably insulation resistance. If, inside a room on fire in 
an wooden historical building, there is not any other element of construction 
with a lower fire resistance, considereing the insulation criteria, then the 
masonry wooden walls will provide some time for fire detection and, eventu- 
ally, fire control before it spreads to other compartments. 


The second observed factor is related with the integrity capacity of the thin 
mortar layer to sustain a permature atack by flames, even when exposed to 
temperatures above 1000ºC. This was in part expected from observations 
carry through following the Chiado fire in 1988, but in this case thicker 
masonry walls were observed and some on then with no wooden section in 
it's interior. 

With a relative thin mortar layer some of the previous beliefs were that such 
a layer would not resist to a temperature exposition higher then 20 minutes 
and after that, all the model will colapse, proving a total usefulness of this 
kind of structure in case of fire. The present conclusions proof that some 
time will be provide for the fire brigade to have arrive before flash-over, even 
if no automatic detection is use, which should not be the case. 


Considering the nature of the buildings in some cities historic centres and, 
above that, the difficulties in security measures improvement due to social 
and economical reasons, the fact that some protection is provide by the exis- 
tant structures, is a good indicator for fire security. 

To conclude is necessary to refer that the historical patrimony is not 
anyone's property but human heritage and therefore should be treated as a 
value to preserve. That is not achieved with simple legislation but should be 
part of a wider education and be reflected in everyday practices. 

The reconstruction process in historic centres considering basic security 
measures for fire prevention is a contribution to that practice and likewise 
should be treated. 


The present work was sponsored by EU - Human Capital and Mobility Pro- 
gramme and the experiments were executed at Cardiff School of Enginee- 
ring, University of Wales. 

The author wishes to acknowledge the contribution of the technical staff of 
the College and of my supervisor, Dr. Anthony Griffiths for his comprehen- 
sion. 

To the colleges in Mouraria Local Committee, and especially to its Co-ordi- 
nator, Argº. Paulo Silva, the author wishes to express his particular recogni- 
tion. 
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Picture 1 Wooden frame in a traditional masonry wall 


Picture 2? Wooden walls under reconstruction in Lisbon historic centre 
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ABSTRACT 


The Aitken's Aº process trans- 
forms a geometric convergent 
sequence into its limit. Taking 
advantage of this fact, some 
features of the convergence 
acceleration processes are 
ilustrated on a conveniently 
chosen example. 

A geometric slowly convergent 
sequence is obtained from the 
application of the Newton- 
Raphson method on the com- 
putation of the root of 
(x-a)'=0,8>1/2. The exact 


root can be computed by the 
Aitken's Aº process using only 
the first three terms of this 
sequence. 

We deliberately omit any argu- 
ment from the iterative proces- 
ses theory. 


RESUMO 


O método A” de Aitken trans- 
forma uma sucessão geomé- 
trica convergente no seu limite. 
Aproveita-se esse facto para 
ilustrar alguns aspectos dos 
métodos de aceleração de con- 
vergência por meio de um 
exemplo convenientemente 
escolhido. 

Aplicando o método de New- 
ton-Raphson ao cálculo da raiz 
de (x-a) =0,6>1/2 obtém- 


-Se uma sucessão geométrica 
lenta. O método A” de Aitken 
permite-nos o cálculo exacto da 
raiz usando apenas os três pri- 
meiros termos desta sucessão. 
Os argumentos próprios da teo- 
ria dos métodos iterativos são 
deliberadamente ignorados. 


PARA QUE SERVEM OS. 
MÉTODOS DE ACELERAÇÃO 
DE CONVERGÊNCIA? 


Mário M. Graça 


Departamento de Matemática. Instituto Superior Técnico. Lisboa 
Original submetido em: 9 de Julho de 1996 


1. INTRODUÇÃO 


Temos sido várias vezes confrontados com a pergunta que serve de título a 
este trabalho. Procuraremos dar a seguir uma resposta aproximada, sob a 
forma dum exemplo de aplicação muito simples, mostrando que há situa- 
ções em Análise Numérica para as quais os métodos de aceleração de con- 
vergência são não apenas desejáveis mas, sobretudo, os processos numéri- 
cos mais adequados à resolução de certos problemas. 


Como ilustração das potencialidades e do interesse dos métodos de acele- 
ração de convergência escolheu-se um método iterativo, o celebérrimo 
método iterativo de Newton-Raphson para a resolução numérica de equa- 
ções reais /(x)=0, no caso em que a função f possui raízes múltiplas. 
Sabe-se que para o cálculo de raizes múltiplas o método iterativo referido, 
se convergente, produz geralmente uma sucessão de valores (x,) conver- 


gindo lentamente para a solução desejada. Os textos de introdução à Aná- 
lise Numérica costumam aconselhar, para este caso, a utilização de outros 
métodos iterativos como por exemplo certas variantes do mesmo método de 
Newton. Todavia a teoria envolvendo a análise de convergência destes 
métodos é, muitas vezes, complexa. 


Neste trabalho, discutimos um exemplo particular para o qual a superiori- 
dade do método de aceleração usado —- o método A? de Aitken [1] — é 
esmagadora. No entanto as ideias aqui utilizadas podem ser ampliadas para 
exemplos mais gerais (o exemplo numérico que tratamos é particular no 
sentido em que conhecemos de antemão a raiz procurada), como seja o 
cálculo numérico de raízes múltiplas de polinómios. O método A* produzirá, 
neste caso, resultados geralmente não tão espectaculares como o que apre- 
sentamos mas, em todo o caso, outras técnicas de aceleração de conver- 
gência existem, com resultados geralmente mais satisfatórios do que os 
obtidos mediante os referidos métodos iterativos. 


Do ponto de vista da aceleração de convergência duma certa sucessão (x,) 


obtida por um qualquer processo iterativo, o método que conduz ao cálculo 
dos seus termos é relativamente indiferente. Se essa sucessão convergir 
para a solução procurada xr, um método efectivo de aceleração de conver- 
gência pode perfeitamente superiorizar-se em relação a todo e qualquer 
método iterativo que escolhamos para a determinação da solução numérica 
da equação f(x)=0. 


No problema de que nos ocupamos na secção 3, o qual consiste na deter- 
minação numérica da raiz da função f(x)=(x-a)”. B>1/2, acontece que a 


sucessão (*,)->« obtida pela fórmula de Newton possui convergência geo- 
métrica lenta (os conceitos e definições necessários à compreensão do 


53 


54 


texto são dados na secção 2). Neste caso particular o método A*, inventado 
por A. C. Aitken em 1926, permite o seu cálculo exacto. 


Como aplicação, num exemplo numérico que tratamos na secção 3, pre- 
tende-se calcular a raíz (de multiplicidade 10) x=1 do polinômio 


plx)= x" -103º +45" —120x +210xº — 252xº + 
-210xº = 1207" +45%º 1044] (1) 


=(x— 1). 


Como se disse, essa raiz pode ser obtida exactamente pelo método A? de 
Aitken. Para o seu cálculo apenas são necessários os três primeiros termos 
da sucessão («| obtidos pelo método de Newton-Raphson. A razão deste 


resultado dramático é, todavia, bem simples. Como se disse, o método A* 
de Aitken permite obter o limite exacto de toda e qualquer sucessão de con- 
vergência geométrica. 


Na secção 2 mostramos também que, para a função f acima, quanto maior 
for o grau de multiplicidade da raiz a calcular, mais lenta é a respectiva 
sucessão de iteradas obtidas mediante aplicação do método de Newton- 
-Raphson. Confirmando o bem conhecido facto de que para polinómios com 
raízes múltiplas, este método pode ser completamente impraticável 
enquanto método iterativo. Note-se porém que do ponto de vista da acelera- 
ção de convergência pouco interessa a maior ou menor lentidão de conver- 
gência da sucessão calculada. No caso presente ela goza duma certa pro- 
priedade (possui convergência geométrica) e isso permite afirmar que 
somos capazes de calcular o seu limite à custa exclusivamente de três 
quaisquer dos seus termos, mediante uma única aplicação da fórmula A* de 
Aitken que relembramos adiante em (5). 

Não é demais sublinhar que, duma forma geral, a análise do tipo de conver- 
gência duma sucessão («| resultante da aplicação dum processo iterativo 
— divergência, convergência linear, convergência geométrica, convergência 
logaritmica, etc. — pode ser obtida exclusivamente a partir dos valores calcula- 
dos (admitindo que essa sucessão goza de propriedades caracteristicas e que 
os seus termos podem ser calculados em número suficiente) abstraindo com- 
pletamente, como já se disse, do método iterativo gerador dessa sucessão. 


No caso presente, inclusivamente, o problema geralmente delicado que con- 
siste na determinação das condições segundo as quais o método de New- 
ton converge para a raiz, qualquer que seja o valor aproximado inicial x, 


não se coloca. Com efeito, mostramos que a sucessão de iteradas (x,) 


converge qualquer que seja o valor inicial +,, sem usarmos nenhuma das 
hipóteses sobre a função f que são dadas habitualmente como condições 
suficientes de convergência do método de Newton. 

Somos assim levados a concluir que os “piores” métodos iterativos de cál- 
culo das raizes duma equação [(x)=0), i.e., os métodos que conduzem a 
sucessões de convergência lenta, muito lenta, ou mesmo muitíssimo lenta, 
podem produzir as “melhores” sucessões possíveis do ponto de vista da 
aceleração de convergência. 


Evidentemente que podemos recorrer ainda a um entre muitos dos outros 
métodos conhecidos de aceleração de convergência [3], (com os mesmos 
resultados no que toca ao problema (1)), no entanto preferimos o referido 
método 4º não só pela sua simplicidade como pelo facto de ele ter come- 
cado recentemente a ser incluido, com frequência, em vários textos universi- 
tários de introdução à Análise Numérica, cerca de 70 anos depois de ter 
sido inventado por A. C. Aitken. 


O autor espera contribuir, através deste modesto trabalho, no sentido de 
fazer notar entre nós o enorme interesse de que se reveste a área de acele- 
ração numérica da convergência de sucessões (que se estende a espaços 
de sucessões escalares, vectoriais, matriciais, etc.). Uma discussão deta- 
lhada sobre os vários métodos de aceleração de convergência, assim como 
descrição das diferentes áreas nas quais eles têm sido aplicados com êxito, 
pode ver-se em [3]. 


2. CONCEITOS BÁSICOS E DEFINIÇÕES 


Nesta secção apresentamos definições e resultados bem conhecidos neces- 
sários à justificação posterior da excelência do método A* quando aplicado 
ao problema particular escolhido. 


* Seja (s,) uma sucessão de valores reais convergente para s. Denota-se 
por (e ja sucessão dos respectivos erros, i.e., (e, )=(s,- s) Diz-se que 


(s,) é uma sucessão de convergência linear se existe o limite 


€, + 


[= lim A (nº) e 


e, 
(a convergência dir-se-á logarítmica quando /=1). Em particular, afirma- 
-Se QUE (s) possui convergência geométrica de razão à, quando a razão 
entre erros consecutivos é constante, i.e. existe Ac R tal que 


Es | 


É (3) 
Por exemplo, as sucessões definidas por s =4+0.999", n=0,1...., 
ucR eu, =1/2'. n=0,1.... são ambas de convergência geométrica. 


Podemos afirmar que a primeira é “mais lenta” que a segunda atendendo 
a que a respectiva razão, calculada mediante (3), é superior à que se 
obtém para a sucessão (u,). 


Uma formulação rigorosa sobre o conceito de “velocidade” de convergên- 
cia duma sucessão pode ver-se em [2]. 


e Dada uma sucessão real (u,) o operador A define-se do seguinte modo: 
Au =u—u, n=0,.. 
Su, <A) k=U3c 


e O teorema seguinte reveste-se de uma importância particular. Ele mostra 
que relativamente a sucessões de convergência linear, se existir o limite 
do quociente dos erros, ele é igual ao limite do quociente das respecti- 
vas diferenças. 


(A notação ( 5) significa lim s, = | Quando ns +). 
Teorema 1 ([6], p. 6) 


Considere-se uma qualquer sucessão real (s,)-»s. Então, supondo 
|+0,1, tem-se: 


—- Sl se e só se 


di (4) 
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O método A? de Aitken é uma transformação de sucessão a sucessão. 
Admitamos que (s ) é uma sucessão para a qual As, =s,-2s,,+ 
+58, £0,VYn, A referida transformação produz uma nova sucessão (1,) 


dada por 


t=s,-————=s, + n=0,1..: (5) 


(a expressão mais à direita em (5) está na origem da designação A* atri- 
buida ao método). 

O teorema seguinte mostra que uma sucessão linearmente convergente 
é transformada pelo processo Aº de Aitken numa sucessão com O 
mesmo limite, e que acelera a convergência da sucessão de partida. 


Teorema 2 ([5], p. 73) 
Seja (s,)-»s uma sucessão de convergência linear tal que As +0, Vn. 


A sucessão obtida pelo método (5) acelera a convergência de (s,), no 


sentido em que 


0) (6) 


Em particular, pode verificar-se por substituição na fórmula (5) de Aitken, 
que toda e qualquer sucessão de convergência geométrica 


s,=b+A, n=01... com beR elâj<1 (para a qual necessariamente 
vn 4's,*0)) é transformada na sucessão constante (1,)=(b). Quer dizer 
que o limite s=» da sucessão (s,) pode calcular-se, mediante a fórmula 


(5), conhecendo apenas os seus três primeiros termos. Obtém-se 1,=5. 


O método de Newton-Raphson, também conhecido pela designação 
abreviada de método de Newton, é um método iterativo que é usado no 
cálculo numérico duma raiz da equação /(x)=0 dada uma sua aproxi- 
mação «,. Às hipóteses mais usuais são que f é uma função continua 
num intervalo fechado onde se encontra a raíz (Única) que se pretende 
calcular. Admite-se ainda que f possui primeira e segundas derivadas 
contínuas no mesmo intervalo ([4], p. 48). O método produz uma suces- 
são (x,) cujos termos são obtidos mediante a fórmula 


dado «,, calcular 


fla) 


ESOM Dado 9 E =0,1.... 7 
Note-se que na fórmula (7) +“. pode sempre calcular-se desde que 
Fix, )+0 vk. Neste caso podemos encarar (7) apenas como uma fórmula 


de recorrência, sem cuidar de saber se se verificam ou não as hipóteses 
habituais de aplicação do método de Newton. Neste sentido, como se referiu 


na introdução, não é o método de Newton em si que nos interessa, mas 


antes uma sucessão particular obtida mediante aplicação da fórmula de 


recorrência (7) iniciada a partir dum certo valor «,. 


3. EXEMPLO DE APLICAÇÃO 


Escolhemos um exemplo para o qual é dispensável toda a teoria envolvendo 
o método de Newton. A análise da convergência da fórmula de recorrência 
(7) é feita directamente sobre a sucessão calculada (x,). 


Considere-se a seguinte função real f definida em RR, onde aceR, 
7 | 
flx ee É E im | 
flxj=(x-a B 3 B (8) 


Para a função f dada em (8), e para uma aproximação inicial x, 5a, a fór- 
mula de Newton (7) escreve-se 


(9) 


Proposição 1 Qualquer que seja a aproximação x, da raíz da equação 
f(x)=0 onde f é a função definida em (8), a fórmula de recorrência (9) pro- 
duz uma sucessão (x,) que converge para a raiz de f (x)=0. 


A sucessão (x,) possui convergência geométrica. A sua convergência é 
tanto mais lenta quanto maior for o valor do parâmetro 5 em (8). 


Demonstração 


Admitamos que lim x, =x (k -» co). Então, por passagem ao limite, de (9) 
resulta necessariamente x=a. Mostremos que, de facto, a sucessão (x,) 


converge. 
Sejam 


=|-—, r=— (10) 
Usando indução matemática sobre k, obtém-se facilmente da expressão (9), 


kl 
mer + k=01... 


Atendendo às hipóteses feitas em (8) com respeito ao parâmetro p, tem-se 
que -1< | < 1 donde 


lima, =rlim SU = E -a (ko) 


Mostremos agora que a sucessão (x,)-»a possui convergência geométrica. 
Aplicando o operador A a ambos os membros de (9) vem 


Ax | 
vk>0 —&=|=|1-=, <=Jel<) 
Ar, B (11) 
e atendendo ao teorema 1, tem-se 
ia Ea? ai À (12) 


x, = 


Como |l|<1 resulta que (x,) é uma sucessão geométrica convergente para 
a. De (11) e (12) concluímos que a sua convergência é tanto mais lenta 
quanto maior for p. 


Como aplicação, passemos ao cálculo da raiz do polinómio dado em (1). 
Tem-se a=1, B=10, 1=9/10. 
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Devido à proposição anterior, possuímos total liberdade de escolha da aproxi- 
mação inicial da raíz do polinómio. Seja x,=-IUO0O. A fórmula de recorrência 
(9) escreve-se 


Ela = (9a, +1) k=01,... 


donde x,=-899.9, x,=-809.81. Apesar destas três aproximações da raiz do poli- 
nómio não sugerirem qualquer indicação sobre a sua grandeza, e a fortiori do 
seu valor exacto, estamos certos de podermos calculá-la exactamente usando 
apenas esses três termos da sucessão (x,) e, com eles, aplicarmos o método 
A? de Aitken (5). Obtém-se, | 

(Ax,) -899.9+ 1000 


n=t————— = 000 —- — nm =, 
EX = 1000 + 2 x 899.9 — 809.8] 
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RESUMO 


O presente artigo analisa os 
caudais de ar de duas salas de 
aula da escola C+S do Crato, 
ventiladas naturalmente com 
recurso a colectores solares a 
ar integrados nas fachadas Sul 
da referida escola. Para levar a 
cabo este estudo foram efectu- 
adas simulações numéricas, 
com recurso ao programa ESP- 
r, para dois meses típicos de 
Inverno e Verão. 

Os resultados das simulações 
enquadram-se com os estudos 
experimentais em curso efectu- 
ados no local, obtendo-se para 
a solução construtiva existente, 
valores compreendidos entre 
120 e 150 mº/h para os cau- 
dais de ar médios no colector. 
Para a satisfação das exigên- 
cias mínimas de qualidade do 
ar sugere-se a utilização de 
dois colectores por sala de aula 
de 2.95mx1.00mx0.20m (em 
lugar de apenas um - e de 
dimensões mais reduzidas - 
como actualmente se verifica). 


ABSTRACT 


The present paper studies the 
airflow performances of two 
classrooms in Crato school 
building, ventilated by natural 
means with resort to solar air- 
collectors integrated in South 
façades. To carry out this 
analysis here numerical simula- 
tions is performed, resorting 
ESP-r program, for two typical 
months in Winter and Summer. 
The results of the simulation fit 
with the experimental survey 
carried out in the local school 
and showed for the existent 
constructive solution, medium 
values of airflow in the solar 


PREDICTED AIRFLOW FROM 
SOLAR AIR-COLLECTORS AT 
CRATO SCHOOL BUILDING 
CASE STUDY 


Paulo R. Pinto (*), A. Canha da Piedade (**) 
Onginal submetido em: 31 de Julho de 1996 


1. INTRODUCTION 


Both indoor air quality and the energy requirements of a building are deeply 
influenced by air infiltration and ventilation. Deficient ventilation may result in 
poor confort, and under extreme conditions, can be injurious to occupants. In 
contrast, unnecessarily high air exchanges rates will present an excessive 
load on heating or cooling systems, resulting the lose of energy or the inca- 
pacity to satisfy thermal requirements. 


The airflow processes are particularly relevant to passive solar buildings, 
since air movement may be the most significant heat transfer mechanism. 
Passive solar heating systems, such as solar air-collectors, may utilise air 
movement mainly by mixed convection to effect heat transfer to a living 
space (for example the classrooms). Air solar panels have, simultaneously, 
two types of air exchange processes: natural and forced convection. Natural 
convection, which arises due to differences in temperature and hence air 
density, creating a pressure difference across the panel that varies with 
height; forced convection, due to wind driving force (wind pressure). 


The Crato school building lies at a latitude of 39.3ºN, a longitude of 7.55ºE, 
and is 300m above mean sea level in the town of Crato. Itis in a rural area 
near Portalegre in a region characterised by significant thermal amplitudes 
(Winter and Summer design temperatures of 0ºC and 34ºC respectively). 
Despite the cool weather during the five months of the year, November to 
March, the region receives substantial sunshine during the winter (total 
amount of 2750 hours/year). The average hemispheric solar radiation inci- 
dent on a vertical south surface for these five months is 490 kWh/m?year; 
however, the annual heating degree-days at 15ºC base temperature are 
908ºC. The weather data show the need for winter heating and a great 
potential for passive solar heating designs. 


2. BUILDING MODELLING 


In the present analysis, a computing environment for building energy simula- 
tion, ESP-r, has been used to predict infiltration, zone-coupled air flow and 
thermal behaviour of the model. It was specified a leakage network to repre- 
sent the distributed air flow paths between nodes or Zzones. The network is 
time dependent to allow the modelling of occupant interaction (timetable of 
the classrooms during the weekdays). An air flow simulation is performed via 
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ESPsim - wind pressure and variable buoyancy driven flows with combined 
heat transfer. 


Building model of the schoolrooms has a total of eleven zones. The criteria 
of the model were based on the prototype classroom of the school building. 
The classroom AN3, aim of our study, is symmetric ffom AN4 (see Fig.1), 
has a total floor area of 67.7 mº and a volume of 227 m”. The classroom has 
west and south-facing glazing with an area of 22% of the pavement. Exterior 
windows of the classrooms have on the outside white plastic blinds (except 
the windows above the roof of the classrooms) in order to prevent thermal 
losses in the cold season and to minimise solar gains in the summer. In the 
cold season, it is admitted that the external blinds are opened during week- 
days between 8.00 a.m. to 6.00 p.m. (legal time - L.T.). In hot weather, 
external blinds are partially closed and joined with interior light coloured cur- 
tains, in order to avoid pernicious solar gains and diffuse the sun's rays. 


section AB 


Fig.1 Diagrammatic representation of the fluid flow network[1] 


The north-south axis of the building was rotated 19.3º (value measured in the 
project) towards the south-east to enable the classrooms to receive gains an 
hour earlier. Building's entrances are reached through transition spaces, 
which prevents the direct contact between indoor and outdoor ambience. The 
south façades of the building overlook a playground and remain free from any 
obstructions that may impinge upon the collection of solar energy. 


air-collectors between 120 and 
150 mº/h, To satisfy the mini- 
mum necessities of air quality it 
is suggested the use of two 
solar air-collectors per class- 
room with 2.95mx1.00mx0.20m, 
instead of the existing one with 
small dimensions. 


O TT ma OS e TT 
= E 


Flat vertical air panels (see Fig.2), each one with size of 2.95mx0.70mx0.20m, 
rely solely on natural convection to transfer solar heat gain. These compo- 
nents of a passive solar feature are an integral element of a building that, 
without mechanical systems, harnesses solar energy. The air panels are 
controlled from inside the rooms by means of simple mechanical commands 
(either inlet or outlet vent may be closed to prevent any unwanted heat out- 
put). The panels are characterised by a non-operable glazing followed by a 
black steel plate for quick absorption and emission of solar energy. As the 
sun's rays pass through the outer glazing and hit the surface of black steel 
plate, this immediately raises its temperature and in turn heats the air trap- 
ped within the glazing and the steel plate. Thermosyphoning could be indu- 
ced by opening the upper and lower vents that allows continuous flow of hot 
air into the classroom during heating periods. The main goal is to increase 
ventilation rates, pre-heating the admitted air and avoiding excessive thermal 
losses. In warm weather the panels will make cross ventilation, since local 
winds are predominantly northerly, linking up with the suction pressure cau- 
sed by the indoor thermal differential (stack effect). 


Fig. 2 View of the solar air-collector 


The type of occupation is discontinuous, from 8.00 a.m. to 6.00 p.m. (L.T.) 
Monday-to-Friday. Internal gains due to the presence of students (each scho- 
olroom is occupied with average 28 students) and the teacher are calculated 
with a total of 1680 W. This is based on the assumption that the heat gain 
rate from a child is 40% of that from an adult, which equals 140 W[2]. 


3. RESULTS AND DISCUSSION 


3.1. Summer period 


Figure 3 illustrates the cross-ventilation rates of schoolroom AN3 induced by 
the air panel with extreme values between 0.3 and 8.8 air change per hour 
(ACH). The minimum value is in harmony with the weekends when the panel 
is not working. During weekdays, the medium value is 1.8 ACH, less than 
the minimum airflow required for the occupation of the classroom (2.3 ACH) 


The airflow exchange due to the panel is between 357 and 1300 m*/h. 
However, the results of the simulation showed in great range of the period, 
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the inverted work of the designed panel, which means that wind pressure is 
preponderant and contrary the buoyant ventilate flow. 


Mode aula] with all comections For day 1 of month 7 (1385) 
to day 31 oF month 7 (1985) 


—— Mode sula3 entering 
—— Mode aula] leaving 
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Fig. 3 Ventilation rates for the classroom AN3 - July(1] 


3.2. Winter period 


Figures 4 and 5 show, respectively, ventilation rates for schoolroom AN3 and 
the solar air-collector. The contribution of the panel to the classroom ventila- 
tion, considering everything closed, duplicates the air renovation rates of the 
classroom for medium values of 0.8 and 0.9 ACH. However, these values 
are not still enough to supply the minimum needs of airflow required for the 
existing occupation of the classroom. 
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Fig. 4 Ventilation rates for the classroom AN3 - January[1] 


— Mode cs3  entering 
— Mode c33 leaving 


5 


Do A) 


“q o % 144 192 240 ME 3% 3 437 430 SM SK 6 672 10 
Time 4h) 


Fig. 5 Ventilation rates for the air panel - January(1] 


The airflow exchange due to the panel is now between 40 and 520 m*/h. 
The medium airflow is 120 mº/h for freefloating regime and can reach 150 
m/h if the schoolroom is under thermostatic temperature of 18ºC. 


Maximum values in the black sheet surface - corresponding to vertical solar 
radiation intensity greater than 500 W/m* - and in the single glass pane are, 
respectively, 77ºC and 45ºC. Indeed, for those values of global radiation the 
average air temperature in the panel is 5 to 10ºC higher than outdoor air 
temperature. 


Despite its positive thermal behaviour, there were some negative feedbacks in 
operating the solar air panel studied in the present analysis. First, it seems to 
be a nuisant daily routine for many students to open and close the vents every 
day. Second, there were some problems with bugs, since the gap between the 
non-operable glazing and the black sheet has become an attractive habitat. 


4. CONCLUSIONS 


Performance analysis of the Crato solar air-collectors do not satisfy air qua- 
lity requirements with resort only the designed existing solution. However, 
this kind of solution can improve indoor air quality with positive thermal 
balance. To sum up, it is suggested to double the number of panels for each 
schoolroom and increase their width (from 0.70 m to 1.00 m). 
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Um programa compacto, potente e fácil de usar porque 
trabalha em ambiente Windows. 

Um programa seguro, com a garantia que só um lider pode 
oferecer. 


A solução perfeita para a engenharia mais exigente. 
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